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Abstract . 


Subroutine  rSTADET"  analyzes  a  structure  to  find  if 
it  is  stable  and  determinate.  The  structure  must  first  be 
idealized  as  a  rigid  body  supported  by  a  system  of  links. 
The  program  checks  on  the  arrangement  and  number  of  the 
support  links .  Flags  are  set  so  that  this  subroutine  may 
be  used  to  select  whether  to  use  simple  statics  or  strength 
of  materials  methods  in  subsequent  analysis  of  the  problem. 
Testing  of  the  subroutine  is  discussed. 

Two  input-output  programs  are  also  described.  One 
calculates  the  reactions  in  the  equivalent  support  linkages 
from  the  resultant  load;  the  other  assumes  that  the  pro¬ 
perties  of  the  links  are  known  and  calculates  the  allowable 
load . 
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Notation . 


A 

A  matrix  of  coefficients 

A ,  B ,  C  ,  D 

The  coefficients  of  the  equation  of  a 
plane 

A 

A  vector 

AAA 

HX>  Ay’  Az 

Rectangular  components  of  A 

B 

A  column  matrix 

B 

A  vector 

Rectangular  components  of  B 

ll>  §2’  etc. 

Vectors  made  up  of  the  direction  cosines 
of  the  support  links,  link  cosine  vectors 

-i J  - j  »  -k 

Indicates  the  ith,  jth  or  kth  of  Bn,  B  , 
etc .  2 

F 

A  Fortran  floating  point,  numerical 
format 

F 

A  vector  force 

— i 

The  ith  vector  force 

Fx>  Fr  Fz 

Rectangular  components  of  F 

^xi,Fyi,Fzi 

Rectangular  components  of  Fj_ 

I 

A  Fortran  integer,  numerical  format 

ij  j.,  k 

Unit  orthogonal  base  vectors  in  the 
x,  y,  z  directions  respectively 

KPROB 

An  array  used  to  Indicate  the  kind  of 
problem  involved 

1 ,  m ,  n 

Direction  cosines  of  a  line 

M 

Vector  moment 

Mx,  My,  Mz 

NOSOL 

Rectangular  components  of  M 

An  array  used  to  indicate  the  reason 
for  the  lack  of  a  solution 

vi 


P].>  »  etc* 

Position  vectors  giving  the  locations 
of  the  points  of  application  of  the 
links  with  respect  to  the  first  point 
of  application 

PLA 

An  array  used  to  store  the  coefficients 
of  the  equations  of  the  planes  determin 
ing  the  lines  of  action  of  the  links 

r 

The  position  vector  of  the  point  of 
application  of  a  force 

rx>  V 

Rectangular  components  of  r 

u 

Poynting's  vector 

W 

A  weight 

X 

A  column  matrix 

X,  y,  Z 

Rectangular  coordinates 

x0  ,y0  >z0 

Coordinates  of  a  given  point 

xi 

The  scalar  result  of  forming  a  dot 
product 

2i 

The  vector  result  of  forming  a  cross 
product 

a 

The  angle  included  between  two  vectors 
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1.  Introduction. 

Digital  computers  have  been  extensively  applied  to 
the  solution  of  mechanics  problems.  However,  these  ap¬ 
plications  have  been  limited  in  scope.  The  computer  has 
been  used  as  a  "super  slide  rule"  to  solve  a  particular 
type  of  problem  which  was  too  time-consuming  or  too  dif¬ 
ficult  previously.  It  is  the  intention  of  Professor  Gilles 
Cantin  to  attack  the  problem  of  the  design  of  structures 
by  a  digital  computer  in  a  general  rather  than  a  specific 
manner . 

One  of  the  first  routines  needed  for  the  Cantin  Pro¬ 
ject  was  a  routine  to  reject  unstable  structures  and, 
perhaps,  make  some  selection  as  to  the  methods  of  analysis 
to  be  used  in  subsequent  tests.  The  subroutine  ’’STADET" 
analyzes  an  idealized  structure  to  find  if  it  is  stable 
and  determinate.  It  sets  flags  which  can  be  used  to  print 
out  a  description  of  the  kind  of  problem  involved.  When 
the  subroutine  finds  the  first  condition  which  would  pre¬ 
clude  the  solution  of  the  problem  by  simple  statics,  a 
flag  is  set  and  control  returns  to  the  input-output  program. 
If  the  "error"  is  caused  by  a  redundant  support  system, 
other  methods  of  analysis  may  be  available  whereas  there 
is  no  solution  for  an  "unstable"  problem. 

The  definition  of  stability  which  is  considered  here 
is  more  restrictive  than  the  question  of  stable  equili¬ 
brium  usually  encountered  in  basic  statics  courses.  Here 
is  meant,  not  the  criterion  of  energy  change  with  small 
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perturbations,  but  the  geometrical  stability  of  infinitely 

strong  members.  As  expressed  by  Kinney, ^ 

A  stable  structure  will  support  any  conceivable 
system  of  applied  loads,  resisting  these  loads 
elastically  and  immediately  upon  their  applica- 
tion  the  strength  of  all  members  and  the  capacity 
of  all  supports  being  considered  infinite. 

Consider  the  two  links  supporting  a  weight  shown  in 


Fig.  1. 


Fig.  1.  Support  links  before  and  after  displacement. 

In  Fig.  la,  with  the  links  horizontal,  there  is  no 
vertical  component  of  the  tension  in  the  bars  to  support  W. 
Even  with  infinitely  strong  foundations  and  members,  the 
geometrical  arrangement  precludes  a  vertical  component  un¬ 
til  some  downward  displacement  of  W  has  occurred  (Fig.  lb). 
The  methods  of  strength  of  materials  are  required  to  find, 
the  equilibrium  position  of  the  structure.  Since  W  is  not 
resisted  "immediately  upon  its  application,"  this  arrange¬ 
ment  fails  the  test  for  stability. 

Both  a  structure  and  its  support  system  must  be 
stable.  The  problem  which  is  examined  here  is  limited 
to  that  of  the  external  stability  of  a  structure.  The 

^J.  Sterling  Kinney,  Indeterminate  Structural  Analysis, 

p.  20,  1957. 
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structure  itself  is  considered  to  be  a  rigid  body  so  that 
only  the  system  of  supporting  links  is  of  concern.  No 
loads  are  imposed  during  the  stability  analysis  since 
equilibrium  is  not  involved. 

However,  in  some  cases  forces  such  as  gravity  must  be 
considered  in  order  to  have  a  structure  which  meets  this 
restricted  definition  of  stability.  This  is  a  familiar 
assumption  in  basic  statics  where,  for  example,  the  gravity 
force  provides  a  reaction  against  upward  motion  of  a  roller 
nest.  It  is  up  to  the  user  of  this  program  to  decide 
whether  a  "gravitational”  restraint  is  sufficient  or  whether 
a  positive  restraint  is  needed. 

The  user  must  first  idealize  his  structure  so  that  the 
subroutine  can  operate  on  it.  Appendix  A,  which  shows  the 
linkage  diagrams  of  some  common  supports,  may  be  helpful  in 
this  regard.  However,  it  is  at  this  point  that  engineering 
judgment  is  required.  Needless  to  say,  if  the  proper  equi¬ 
valent  support  system  is  not  employed,  no  reliance  can  be 
placed  on  any  answers  which  are  produced. 

Fig.  2  shows  the  decisions  which  are  made  by  the  sub¬ 
routine.  The  titles  within  the  blocks  indicate  the  basic 
sections  of  the  subroutine.  These  titles  are  roughly  indica¬ 
tive  of  the  purpose  of  the  sections  though  more  information 
is  usually  produced  than  the  title  would  lead  one  to  believe. 

For  input  the  subroutine  requires  the  direction  cosines 
of  the  support  linkages,  the  coordinates  of  the  points  of 
application  of  the  links  and  the  number  of  links.  Appendix  H 
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Figo  2  Simplified  Flow  Chart 
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shows  how  the  input  must  be  arranged  to  be  read  in  by  the 
input-output  programs  of  Appendix  C.  The  output  consists 
of  the  two  flag  arrays,  KPROB  and  NOSOL. 

The  subroutine  first  considers  the  points  of  appli¬ 
cation  of  the  supporting  links.  If  these  points  are  along 
a  line,  the  question  of  whether  it  is  a  two-  or  three- 
dimensional  problem  depends  on  the  orientation  of  the  links. 
If  the  points  are  in  a  plane,  the  links  must  all  be  in  that 
same  plane  for  a  two-dimensional  problem  to  be  involved.  If 
the  points  are  randomly  located  in  space,  only  a  three- 
dimensional  problem  can  be  involved.  The  two-dimensional 
problem  requires  three  linkages  for  stability;  the  three- 
dimensional  problem,  six.  If  there  are  fewer  linkages,  an 
"unstable"  error  print-out  is  made;  if  more,  an  "indeter¬ 
minate"  one  is  made. 

When  it  has  been  determined  that  the  correct  number  of 
support  linkages  is  involved,  the  orientation  of  the  links 
is  considered.  The  following  restrictions  apply  to  the 
three-dimensional  case: 

1.  All  the  links  may  not  be  parallel. 

2.  They  may  not  lie  in  parallel  planes. 

3.  The  axes  of  the  links  may  not  intersect  one 
straight  line  since  limited  rotation  would 
be  possible  about  that  line. 

4.  The  axes  of  the  links  may  not  intersect  in  a 
point  or  limited  rotation  is  possible  about 
that  point. 

5.  Four  of  the  six  links  may  not  be  parallel  or 
intersect  in  any  one  point  for,  in  these  cases, 
a  straight  line  can  be  found  through  the  point 
(or  infinity)  which  intersects  the  axes  of  the 

two  remaining  links  and  restriction  3*  is  violated. 

In  the  two-dimensional  case  only  restrictions  1.  and 


4 .  apply . 
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Where  the  links  are  all  parallel  to  a  given  line, 
limited  motion  is  possible  in  a  direction  perpendicular 
to  this  line. 

In  Fig.  3?  where  the  pairs  of  links  are  in  parallel 


Fig.  3-  Arrangement  with  links  in  parallel  planes. 

planes,  some  motion  is  possible  in  a  direction  perpen¬ 
dicular  to  these  planes  before  a  restoring  component  is 
developed.  In  Fig.  4  some  rotation  about  the  intersection 


* 

Fig.  4*  Arrangement  with  intersecting  lines  of  action 
of  links. 

of  the  lines  of  action  of  the  links  is  necessary  before 
a  restoring  couple  is  developed. 

The  following  section  and  the  appendices  describe 
the  "STADET’1  subroutine,  its  auxiliary  subroutines  and 
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the  Input-Output  Programs  which  were  used  to  test  "STADET" . 
Fortran  language  was  used  for  all  the  programs.  The  testing 
was  done  on  a  Control  Data  Corporation  1604  digital  computer. 
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2.  Method. 

Subtitles  indicate  the  sections  of  the  flow  charts 
(Figs.  5  through  9)  and  the  "STADET''  list  (Appendix  J) 
under  discussion.  Circled  numbers  on  the  flow  charts  are 
used  for  two  purposes.  The  first  is  to  facilitate  the 
comparison  with  the  list  of  Appendix  J;  the  second  is  to 
indicate  a  transfer  to  another  section  of  the  program 
when  the  use  of  a  line  to  show  the  path  would  be  cumber¬ 
some  . 

The  subroutine  only  operates  up  to  the  detection  of 
the  first  error.  The  appropriate  flag  is  then  set  and  a 
return  is  made  to  the  Input-Output  Program. 

a.  C0PLANAR1 

This  section  of  the  subroutine  generates  the 
position  vectors  which  give  the  locations  of  the  points 
of  application  of  the  links.  It  checks  to  see  whether 
these  points  are  coincident,  lie  on  a  line,  are  in  a  plane 
or  are  randomly  located  in  space.  In  the  case  where  the 
points  of  application  are  in  a  plane,  this  section  of  the 
program  checks  to  see  if  the  links  are  in  the  same  plane 
and  so  determines  whether  a  planar  problem  is  involved. 

After  deciding  whether  a  two-  or  three-dimensional 
problem  is  involved,  this  section  counts  the  number  of  links 
to  find  out  if  the  number  is  sufficient  for  stability. 

The  computer  operates  with  a  binary  approximation  to 
a  decimal  number  when  it  is  used  in  the  floating  point- 
mode.  Each  calculation  is  approximate  and  round-off  errors 


Fig.  5  Flow  Chart  of  C0PLANAR1 


can  accumulate.  For  this  reason  the  "zero’’  which  must  be 
used  in  the  program  is  not  the  mathematical  zero  but  an 
allowable  tolerance.  The  first  portion  of  this  section 
examines  the  coordinates  of  the  points  of  application  and 
determines  the  value  of  ’’zero”  to  be  used  throughout 
"STADET".  Appendix  L  discusses  the  method. 

Vectors  are  generated  between  the  first  point  of 
application  and  the  other  points.  These  vectors  are  ex¬ 
amined  to  find  the  first  non-zero  one.  If  none  of  them 
is  non- zero,  there  is  only  one  point  of  support,  the  struc¬ 
ture  is  unstable,  and  a  return  is  made  to  the  Input-Output 
Program. 

If  a  non- zero  vector  is  found,  the  subroutine  deter¬ 
mines  whether  it  is  the  last  one  or  not.  If  it  is  not 
the  last,  the  other  vectors  are  examined  to  find  whether 
any  others  are  non-zero.  If  the  first  non-zero  vector  is 
the  last  of  the  array  of  position  vectors,  or  if  all  after 
the  first  are  zero,  the  points  of  application  of  the  links 
lie  on  one  line. 

Having  determined  that  the  points  of  application  are 
on  a  line,  the  subroutine  considers  the  number  of  links. 

A  three-dimensional  linkage  system  is  unstable  because  the 
links  intersect  a  line  and  only  a  planar  problem  has  a 
unique  solution.  Therefore,  if  there  are  more  than  three 
links,  the  system  is  redundant.  If  there  are  less  than 
three  links,  the  system  is  unstable.  If  the  number  of 
links  is  not  satisfactory,  the  appropriate  flag  is  set. 
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Where  a  second  non- zero  vector,  Pj ,  is  found,  its 
cross  product  with  the  first,  P-^,  is  formed.  The  cross 
product  vector,  is  examined  to  find  out  if  it  is  zero. 

If  it  is  zero,  P^_  and  Pg  are  parallel.  Since  Pj_  and  P£ 
have  the  same  initial  point,  this  mades  them  coincident 
and  indicates  that  the  first  three  points  of  application 
are  on  a  line.  The  subroutine  then  checks  to  see  if  all 
of  the  position  vectors  have  been  examined.  If  they  have 
not  all  been  examined,  the  search  for  another  non-zero 
vector  continues.  If  they  all  have  been  examined  and  all 
of  the  points  of  application  are  found  to  lie  on  one  line, 
an  investigation  of  the  number  of  supporting  links  is  made, 
as  above. 

If  a  cross  product  vector,  as  found  above,  Is  not 
zero,  a  check  is  made  of  the  number  of  supporting  links. 

If  there  are  only  three  links,  there  are  only  two  position 
vectors  which  (since  they  are  not  parallel)  determine  a 
plane.  If  there  are  more  than  three  links,  there  are  more 
than  the  two  vectors  that  were  crossed  to  form  Xq • ~  To 
check  whether  the  points  of  application  are  in  a  plane, 
it  is  necessary  to  dot  these  other  vectors  with  Xq*  The 

dot  product,  X^ ,  a  scalar,  is  compared  with  zero.  If  all 
the  dot  products  are  zero,  the  points  of  application  lie 
in  one  plane. 

With  the  points  of  application  of  the  links  in  a  plane 
a  stable  system  is  possible  with  either  a  two-  or  three- 
dimensional  linkage  system.  If  all  the  links  lie  in  the 
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plane  that  contains  the  points  of  application  of  the  links, 
the  problem  is  one  in  two  dimensions.  To  check  for  a  two- 
dimensional  problem,  the  vectors  made  up  of  the  direction 
cosines  of  the  links  are  dotted  with  and  the  dot  products 
are  compared  with  zero.  Since  the  link  cosine  vectors  have 
one  point  in  the  point-of-application  plane,  a  zero  dot 
product  indicates  that  the  vector  lies  entirely  in  that 
plane.  If  a  planar  problem  is  what  is  involved,  i.e.,  all 
the  dot  products  are  zero,  the  proper  number  of  links  is 
checked  for,  as  before. 

When  the  points  of  application  are  not  in  a  plane  or 
when  the  links  are  not  in  the  point-of-application  plane, 
a  three-dimensional  problem  is  involved  and  six  links  are 
required  for  stability.  Suitable  error  flags  are  set  if 
there  are  not  six  supports. 

At  this  point  it  has  been  determined  whether  a  two- 
or  three-dimensional  problem  is  involved  and  that  the 
correct  number  of  supporting  links  is  involved  for  the 
type  of  problem.  Next,  instability  due  to  incorrect 
orientation  of  the  links  is  investigated. 

b .  PARALLEL! 

This  section  checks  to  see  if  the  supporting 
links  are  parallel  or  lie  in  parallel  planes.  The  first 
link  cosine  vector  is  crossed  with  each  of  the  succeeding 
ones  and  the  resulting  cross  product  vectors  are  compared 
with  zero.  If  all  of  the  cross  products  are  zero,  all  of 
the  links  are  parallel. 
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Fig°  6  Flow  Chart  for  PAPALLEL1 
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The  second  portion  of  this  section  is  concerned  with 
discovering  whether  a  situation  similar  to  the  one  shown 
in  Fig.  3  exists  or  not,.  In  the  above  test  a  non- zero 
cross  product  will  result  if  two  of  the  links  are  not 
parallel o  (If  the  situation  is  like  that  of  Fig.  3,  the 
cross  product  vector  is  in  the  direction  of  motion  of  the 
rigid  body.)  If  all  of  the  links  are  perpendicular  to 
this  cross  product  vector,  the  case  is  indeed  similar  to 
that  of  Fig.  3»  This  section  checks  whether  the  links 
are  perpendicular  by  dotting  each  of  the  link  cosine 
vectors  with  the  cross  product  vectors  and  comparing  the 
resultant  scalar  with  zero.  If  all  of  the  dot  products 
are  zero,  all  of  the  link  cosine  vectors  are  in  parallel 
planes.  If  all  are  in  parallel  planes,  the  appropriate 
flag  is  set;  if  not,  the  program  goes  on  to  the  next 
section . 

c .  PARALLEL2. 

Even  though  all  of  the  links  are  not  parallel, 
the  supporting  linkage  system  is  unstable,  by  criterion 
3.,  if  four  of  the  six  links  are  parallel.  Therefore, 
a  separate  check  is  made  for  this  possibility.  The  cross 
product  test  is  used  along  with  a  counting  system.  If  the 
linkage  system  passes  this  test,  the  remaining  check  for  a 
three-dimensional  system  is  done  by  C0NCURRENT1. 

d.  PLANAR1 

This  section  of  the  subroutine  checks  the  three- 
link  cases  to  see  if  the  links  are  in  a  plane  and  if  the 
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Figo  7  Flow  Chart  for  PARALLELS 
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Figo  8  Flow  Chart  of  PLANAR1 


-^return] 

s  parallel 
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links  are  parallel.  In  Fig.  2  there  are  two  paths  from 
C0PLANAR1  to  PLANAR!.  The  first  of  these  three-link  cases 
has  the  points  of  application  on  a  line.  In  the  other 
case  the  points  of  application  are  in  a  plane  and  it  is 
known  that  the  links  are  in  the  same  plane.  For  this 
second  case  there  is  some  duplication  of  testing,  but 
the  extra  operating  time  is  negligible  and  there  was  some 
convenience  in  programming  this  way. 

This  section  crosses  the  first  link  cosine  vector 
with  the  second.  If  the  cross  product  is  zero,  the  first 
is  crossed  with  the  third  link  cosine  vector.  If  this 
cross  product  is  also  zero,  the  three  links  are  parallel 
and  the  structure  is  unstable.  If  the  second  cross  product 
is  not  zero,  it  is  dotted  with  the  second  link  cosine 
vector.  If  the  dot  product  is  zero,  all  three  link  cosine 
vectors  are  in  the  same  plane,  a  stable  planar  problem  is 
possible  and  the  program  goes  to  C0NCURRENT1.  If  the  dot 
product  is  not  zero,  the  three  links  are  not  coplanar  and 
the  system  is  unstable. 

When  the  first  cross  product  is  not  zero,  it  is  dotted 
with  the  third  link  cosine  vector  and  the  dot  product  is 
compared  with  zero.  If  the  dot  product  is  not  zero,  the 
system  is  unstable,  as  above.  If  the  dot  product  is  zero, 
the  program  goes  to  CONCURRENT!. 

e.  C0NCURRENT1 

The  remaining  check  is  to  find  if  the  lines  of 
action  of  the  links  intersect  at  a  point.  In  order  to 
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Figo  9  Flow  Chart  of  CONCURRENT! 


be  able  to  operate  mathematically,  it  is  necessary  to 
convert  the  link  cosine  vectors  to  lines  in  space.  The 
line  form  which  was  chosen  was  that  of  the  intersection 
of  two  planes.  (See  Appendix  G  for  discussion  of  sub¬ 
routine  PLANE.) 

All  of  the  planes  are  determined  and  the  coefficients 
are  stored  in  the  array  PLA.  The  equations  in  PLA  are 
examined  four  at  a  time  but  with  an  overlap,  i.e.,  equations 
one  through  four  are  examined  first,  then  equations  three 
through  six;  and,  if  a  six-link  problem  which  gives  12 
equations  is  involved,  the  process  continues  until  ail 
12  have  been  examined. 

The  three-link  system  must  have  at  least  one  inter¬ 
section  or  the  links  would  all  be  parallel,  a  possibility 
which  has  previously  been  eliminated.  In  a  six-link 
problem  at  least  four  of  the  lines  of  action  of  the  links 
must  intersect  to  make  the  linkage  system  unstable.  There¬ 
fore,  the  system  of  examining  pairs  of  lines  whose  equations 
are  adjacent  in  the  array  PLA  must  find  the  invalidating 
intersection  and  a  more  complicated  checking  system  is  not 
required . 

The  coefficients  of  four  equations  are  treated  as 
though  they  were  the  coefficients  of  four  simultaneous 
equations  in  four  unknowns.  They  undergo  the  matrix 
singularity  test  of  REAC1  (See  Appendix  F).  If  the  de¬ 
terminant  of  the  matrix  is  not  zero,  the  set  of  equations 
is  consistent  and  the  first  two  lines  do  not  intersect. 
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If  the  next  overlapping  set  of  four  is  also  consistent, 
three  lines  do  not  intersect  and  the  link  support  system 
is  satisfactory  „  This  would  complete  the  checks  by  " STADET" . 

On  the  other  hand,  if  the  determinant  is  zero,  the  set 
of  equations  is  dependent  or  inconsistent  and  the  two  lines 
intersect  or  are  parallel.  If  the  lines  intersect,  the 
point  of  intersection  can  be  found  by  solving  three  of  the 
equations  simultaneously.  If  the  lines  are  parallel,  the 
four  combinations  of  these  four  equations,  taken  three  at 
a  time,  are  all  inconsistent.  No  account  is  taken  of  this 
parallelism  since  it  is  already  known  that  less  than  four 
links  are  parallel.  The  next  overlapping  set  of  four  equa¬ 
tions  is  then  considered. 

When  an  intersection  is  found,  the  coordinates  are 
substituted  into  the  equations  of  all  the  planes.  The 
number  of  equations  which  are  satisfied  is  counted.  In 
the  case  of  a  planar  problem,  if  six  are  satisfied,  all 
the  lines  of  action  intersect  and  the  linkage  system  is 
unstable.  In  the  three-dimensional  case,  if  eight  or 
more  are  satisfied,  the  linkage  system  is  unstable.  Where 
instability  exists,  the  appropriate  flag  is  set. 

This  completes  the  checks  made  by  "STADET”.  The 
rigid  body  is  supported  by  a  system  of  links  whose  re¬ 
actions  can  be  determined  by  simple  statics  or  the  appro¬ 
priate  flag  has  been  set  in  NOSOL. 

The  one  remaining  possible  error,  that  of  not  having 
the  load  in  the  correct  plane  for  a  planar  problem,  is  not 
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strictly  of  concern  in  deciding  whether  the  system  of 
links  is  stable  and  determinate .  The  check  for  this 
error  is  discussed  in  Appendix  C,  Input-Output  Programs. 
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3.  Testing  of  Program. 

In  order  to  devise  a  logical  system  of  tests  the  sim¬ 
plified  flow  chart  of  Fig.  2  was  used.  This  shows  the 
questions  that  the  various  sections  of  the  subroutine 
answer.  In  all,  there  are  25  distinct  paths  that  can  be 
followed.  Only  four  of  these  paths  lead  to  a  satisfactory 
problem,  i.e.,  one  which  can  be  solved  by  simple  statics. 
Fig.  10  lists  the  types  of  problems  which  were  solved.  The 
complete  print-outs  of  the  solved  problems  are  collected 
in  Appendix  K. 

Problems  26  through  29  are  the  complements  of  the 
satisfactory  problems  of  the  first  25.  The  "bar  reactions 
output  of  TEST1A  was  used  as  the  "allowable  bar  loads"  of 
TEST1B.  The  output  of  TESTIB,  the  "allowable  load  terser* 
should  be,  and  is,  the  "resultant  torsor"  of  TEST1A. 

In  problems  11  and  17,  where  there  are  more  than  six 
links,  the  storage  process  is  such  that  the  last  link  is 
stored  on  top  of  the  first.  This  results  in  a  nonsense 
print-out  of  the  structure  data.  In  problem  11,  for  ex¬ 
ample,  the  first  link  has  two  direction  cosines  of  one. 

In  the  use  of  TEST1A  to  solve  a  planar  problem  it 
is  necessary  to  consider  the  possibility  of  having  an 
invalid  problem  because  the  load  is  not  in  the  same  plane 
as  the  structure.  Since  this  possibility  is  caused  by  an 

''Torsor  --  A  coined  word  (combination  of  torque  and 
tensor)  used  to  indicate  a  generalized  resultant  which 
is  composed  of  a  couple  and  a  force. 
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List  of  Problems 


a  c 

o  o 

p  P 

oj  cti 

o  o 

-H  P 

tH  rH 

CX  cx, 

O-a  0-4 

o5  ctf 

Ch  0  Ps 

GO  O  U 

®HS  05 

CO  rH  CO  £3 

P  P  ctf 

C  05  0  H 

°H  °H  Oa 

G  d  O  O  n 

o<  o  o  Comments 


1 

X 

Points  are  concurrent 

2 

X 

X 

Less  than  three  links  (supports) 

3 

X 

X 

More  than  three  links 

4 

X 

X 

Links  are  parallel 

5 

X 

X 

Links  are  not  in  one  plane 

6 

X 

X 

Lines  of  action  of  links  intersect 

7 

X 

X 

More  than  three  links 

a 

X 

X 

Links  are  parallel 

9 

X 

X 

Lines  of  action  of  links  intersect 

10 

X 

X 

Less  than  six  links 

11 

X 

X 

More  than  six  links 

12 

X 

X 

Links  are  parallel 

13 

X 

X 

Links  are  in  parallel  planes 

14 

X 

X 

Four  of  six  links  are  parallel 

15 

X 

X 

Lines  of  action  of  links  intersect 

16 

X 

Less  than  six  links 

17 

X 

More  than  six  links 

ia 

X 

Links  are  parallel 

19 

X 

Links  are  in  parallel  planes 

20 

X 

Four  of  six  links  are  parallel 

21 

X 

Lines  of  action  of  links  intersect 

22a 

X 

X 

22b 

X 

X 

Resultant  force  not  in  plane  of 
structure 

23a 

X 

X 

23  b 

X 

X 

Resultant  force  not  in  plane  of 
structure 

23  c 

X 

X 

Resultant  force  not  exactly  in 
plane  of  structure 

24 

X 

X 

25 

X 

26 

X 

X 

27 

X 

X 

Oblique  plane  problem 

23 

X 

X 

Stable  spatial  problem 

29 

X 

Stable  spatial  problem 
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error  on  the  part  of  the  one  setting  up  the  problem  and 
is  not  strictly  a  factor  in  the  stability  of  the  structure, 
testing  for  this  possibility  was  made  part  of  the  Input- 
Output  Program.  The  testing  of  a  planar  problem  therefore 
requires  one  step  beyond  T'STADET" .  The  subcases  of  prob¬ 
lems  22  and  23  show  the  results  of  such  tests. 

It  was  not  possible  to  get  an  accurate  estimate  of 
the  time  required  to  solve  a  problem.  It  takes  roughly 
one  second  to  solve  an  unsatisfactory  problem  and  about 
twice  as  long  to  solve  a  satisfactory  one.  This  differ¬ 
ence  is  not  in  the  computation  time  but  in  the  time  re¬ 
quired  to  write  the  output  on  the  magnetic  tape.  In  fact, 
when  a  series  of  problems  is  being  solved,  only  a  slight 
pause  can  be  detected  between  the  end  of  writing  one 
answer  and  the  beginning  of  the  output  for  the  next.  To 
these  computing  and  writing  times  must  be  added  the  time 
to  compile  the  program.  This  is  about  75  seconds  for 
TEST1A  and  about  15  seconds  for  TEST1B. 
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4.  Conclusion o 

When  this  project  was  undertaken,  it  was  understood 
that  this  subroutine  would  eventually  be  used  as  part  of 
a  much  more  ambitious  program .  The  objective  of  providing 
a  routine  which  decides  whether  the  structure  may  be 
analyzed  by  the  methods  of  simple  statics  has  been  achieved. 

For  the  purposes  of  this  routine  a  decision  was  made 
to  have  storage  for  only  six  equivalent  links,  the  maximum 
number  compatible  with  a  solution  obtainable  from  statics. 
This  is  the  only  restriction  which  must  be  removed  in 
future,  expanded  use.  The  arrays  KPROB  (kind  of  problem} 
and  NOSOL  (no  solution)  were  set  up  to  be  used  as  keys  to 
select  the  proper  methods  of  attack.  If  N0S0L{4)  is  set- 
to  one,  a  redundant  structure  is  involved  and  more  power= 
ful  methods  must  be  used  to  effect  a  solution.  KPROB 
indicates  whether  the  simplification  of  a  planar  problem 
is  valid. 

While  if  is  probable  that  a  bona  fide  planar  problem 
would  be  set  up  in  one  of  the  coordinate  planes,  it  might 
be  advisable  to  provide  a  subroutine  to  transform  the 
coordinate  system  so  that  an  oblique  plane  problem  is 
transferred  to  a  coordinate  plane.  It  may  or  may  not  be 
desirable  to  transform  the  results  back  to  the  original 
coordinate  system. 

The  next  step  in  the  generalized  analysis  of  structures 
by  digital  computers  could  be  along  one  of  several  paths. 
Since  this  subroutine  is  concerned  with  only  the  external 
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stability  of  the  structure,  a  subroutine  to  examine  the 
internal  stability  needs  to  be  written „  The  simplifica¬ 
tion  of  considering  only  space  frames  might  yield  a  good 
first  approximation o 

Another  path  would  be  to  either  use  a  code  or  have 
an  examining  type  of  subroutine  to  decide  the  type  of 
problem  and  then  jump  to  one  of  a  set  of  subroutines 
designed  to  solve  a  given  class  of  problems  * 

It  is  felt  that  this  is  one  small  step  in  removing 
some  of  the  drudgery  from  engineering  design. 
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APPENDIX  A 


Linkage  Diagrams  of  Common  Supports 


A  link  is  a  member  capable  of  supporting  only  an 
axial  force.  It  has  spherical  hinges  at  both  ends  so 
that  it  cannot  transmit  any  moments.  In  the  drawings 
below  the  link  is  shown  thus:  ,  where  the  dot  indi¬ 

cates  the  hinge  at  the  point  of  application.  The  other 
end  of  the  link  is  to  be  understood  to  have  a  hinge  which, 
in  turn,  is  attached  to  a  frame  of  reference. 

As  discussed  in  Section  1,  Introduction,  some  of 
these  supports  are  not  the  equivalent  of  a  complete  link. 
In  these  cases  it  is  the  user's  responsibility  to  either 
modify  the  support  or  to  make  sure  that  the  reaction  is 
in  the  allowed  direction. 


Smooth  plane  support 


Because  friction  is  neglected, 
a  linkage  system  which  only 
opposes  vertical  motion  is 
obtained.  The  link  shown  may 
only  have  a  compressive  reaction. 


Flexible  cord 


The  reaction  is  along  the  cord. 
Only  a  tensile  reaction  is 
allowable . 
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Ball  and  socket 


A  ball  and  socket  is  the 
equivalent  of  three  links 
acting  at  a  point. 


Built-in  support 


BODY 


A  built-in  support  is  the 
equivalent  of  six  links. 

The  transparent  body  is 
shown  supported  by  one  of 
many  stable  and  determinate 
arrangements  of  six  links. 
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APPENDIX  B 


Mathematical  Tests 


The  examination  of  the  system  of  supporting  links  is 
accomplished  by  the  application  of  a  few  simple  mathemati¬ 
cal  tests  to  the  vectors  which  are  part  of  or  are  generated 
from  the  input  data. 

A.  The  dot  or  inner  product  of  two  vectors  A  and  B, 
neither  of  which  is  zero,  is  zero  if  and  only  if  A  and  B 
are  perpendicular. 

A  •  B  =  [A|  |B|  cos  a  where  a  =  included  angle 


A  •  B  =  AXBX  +  AyBy  +  AZBZ 


subscripts  indicate  x, 
y,  z  components  of  A 
and  B 


Bo  The  cross  or  outer  product  of  two  vectors  A  and 
B,  neither  of  which  is  zero,  is  zero  if  and  only  if  A  and 
B  are  parallel  or  collinear. 

A  x  B  =  |A|  | B |  sin  a  u  where  u  is  a  unit  vector  per¬ 

pendicular  to  both  A  and 
B  and  with  them  forming 
a  right-handed  set 


A  x  B  = 


i  i  k 

^x  Ay  A z 

B  B  B 
x  y  z 


i,  2  and  k  are  unit,  orth¬ 
ogonal  base  vectors  in  the 
x,  y,  z  directions,  re¬ 
spectively 


C.  The  question  of  whether  of  not  lines  intersect  is 
solved  by  examining  the  equations  of  the  planes  determining 
the  .lines o  If  a  set  of  four  equations  determining  a  pair 
of  lines  is  inconsistent  or  dependent,  the  lines  are  parallel 
or  intersect o  The  "matrix  singularity  test"  of  REAC1  (see 
Appendix  F)  is  used  to  find  if  the  equations  are  dependent 
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or  inconsistent* 

If  the  lines  intersect,  the  point  of  intersection 
can  be  found  by  solving  three  of  the  four  equations  simul¬ 
taneously* 


31 


APPENDIX  C 

Input-Output  Programs 


These  control  programs  provide  for  reading  in  data, 
operating  the  required  subroutines  and  printing  results. 

In  both  cases  the  first  step  is  reading  in  the  data  on 
the  equivalent  linkages.  Then  "STADET"  is  called  and  the 
stability  and  determinateness  of  the  support  system  is 
examined.  The  rest  of  the  data  for  that  problem  is  then 
read  in  regardless  of  the  results  of  "STADET".  The 
intermediate  operations  of  the  two  programs  will  be  con¬ 
sidered  separately. 

TEST1A 

If  the  structure  is  stable  and  determinate,  the  re- 
sultant  of  the  forces  and  couples  acting  on  the  structure 
is  calculated.  This  establishes  the  B  vector  in  the  matrix 
equation,  Ax  =  B.  The  A  matrix,  that  array  which,  when 
used  to  premultiply  the  link  reaction  vector,  x,  gives  B, 
is  found  next.  If  there  are  six  links,  A  is  six  by  six; 
the  six  equations  represented  are  independent  and  the  values 
of  the  link  reactions  are  found  by  subroutine  REAC1. 

In  the  case  of  a  planar  problem,  "STADET''  has  examined 
the  structure  and  found  it  satisfactory,  but  it  remains  to 
be  seen  whether  the  load  is  applied  in  the  correct  plane. 
Therefore,  checks  are  made  to  see  that  the  resultant  force 
vector  is  in  the  plane  of  the  structure  and  that  the  re- 
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sultant  couple  vector  is  perpendicular  to  that  plane . 

If  there  are  only  three  links,  A  is  six  by  three; 
only  three  independent  equations  exist  among  the  six, 
and  it  is  necessary  to  solve  various  sets  of  three  equa¬ 
tions  until  a  consistent  set  is  found  (See  Appendix  F, 
REAC1).  When  a  consistent  set  is  found,  the  results  are 
printed  out  as  discussed  below. 


TEST1B 

If  the  structure  is  stable  and  determinate,  the  A 
matrix  is  calculated.  The  allowable  loads  in  the  various 
links  are  knowns  in  this  problem.  Therefore,  the  x  vector 
is  multiplied  by  A  giving  the  allowable  load  on  the  struc- 
ture.  In  this  case  there  is  no  question  about  the  plane 
of  the  resultant  load  as  in  TEST1A. 

Output 

In  both  programs  the  intermediate  calculations  are 
suppressed  if  the  structure  is  not  stable  and  determinate. 
The  printout  then  consists  of  the  input  structure,  a 
description  of  the  problem  type  and  the  error.  If  the 
structure  is  satisfactory,  the  input  for  the  intermediate 
calculations  and  their  results  are  printed.  Since  the 
case  of  coincident  support  points  would  be  an  unusual 
blunder,  a  special  error  printout  was  not  used.  Since 
"STADET"  is  not  really  operated  in  this  case,  the  print¬ 
outs  are: 
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NOT  A  PLANAR  PROBLEM  and 

STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS 
INTERSECT  A  LINE 

(See  problem  1.) 

Both  programs  reset  conditions  so  as  to  solve  the 
next  problem  and  then  continue  until  all  solutions  have 
been  made.  (See  Appendix  H  for  problem  format.) 
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PROGRAM  TEST1A 

DIMENSION  VB  (  6 , 3)  , B( 6,3) , FORCE ( 100,3) ,CORD( 1 0 C, 3 ) , COUPL ( 50 , 3 ) , 
1KPRCB15) , NOSOL (6) ,RF0R(3) , RCCPL ( 3 ) , BREAC ( 6 ) , CCMX (6 , 7 ) ,PS16,7), 

2 AR ( 3 )  , BR ( 3 ) , CROSS ( 3 )  , VP ( 5 , 3 ) , PO INT ( 3 , 4 ) , PLA  ( 1 2 ,4 ) , PL ( 4  )  , PM ( 4 ) 
COMMON  VB,B, FORCE , CORD, COUPL ,NFO,NCC, NB AR ,KPROB, NO SOL , RFOR , RCOPL , 

1 BRE  AC , COMX » VP 
5  FORMAT (15) 

10  FORMAT! Ii, 12) 

19  FORMAT! 3F10. 8, 3F1  4.7) 

20  FORMAT ( 6F 1 2 . 6 ) 

30  FORMA T13F14.7) 

90  FORMAT!/////) 

400FORMA  T ( 80H 

1  ///) 

1C1  FORMA T ! 5X ,  48HEXA MINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE//) 

1 110F0RMAT! 1 1X,25HDIRECTION  COSINES  OF  BARS,21X, 

1 2 1 HPO I  NTS  OF  APPLICATION) 

121  FORMAT  (8X, 1HL, 12X, 1HM, 12X, 1HN, 13X, 1HX, 16X, 1HY, 16X, 1HZ) 

129  F0RMAT(3!2X,F11.7),3(3X,F14.7)) 

91  FORMAT!//) 

51COFORMAT(5X,  61HSTRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  C 


10NCURRENT) 


VjJ 

O 


5100FORMATI5X,  61HSTRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  C 
1 ONCURRENT ) 

51 10F0RMAT(5X,59HSTRJCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  PA 
1RALLEL) 

512  FORMAT (5Xf50H STRUCTURE  UNSTABLE  INSUFFICIENT  NUMBER  OF  SUPPORTS) 

513  FORMAT (5X»49HSTRJCTURE  INDETERMINATE  EXCESS  NUMBER  OF  SUPPORTS) 
5140F0RMAT (5X»63HSTRJCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERS 

1 ECT  A  LINE) 

515  F0RMAT(5X»27HREACTI0N  MATRIX  IS  SINGULAR) 

516  FORMAT (5X»28HSUPPORT  POINTS  ARE  ON  A  LINE) 

517  FORMAT (5X,29HSUPP0RT  POINTS  ARE  IN  A  PLANE) 

518  FORMAT (  5X  * 14HPLANAR  PROBLEM) 

519  FORMAT! 5 X , 4 3HRES J LT ANT  LOAD  IS  NOT  IN  PLANE  OF  STRUCTURE) 

521  FORMAT (5X»20HN0T  A  PLANAR  PROBLEM) 

522  FORMAT (5X»51H STRUCTURE  UNSTABLE  SUPPORTS  ARE  IN  PARALLEL  PLANES) 
100  FORMA  T ( 5X »  48HFIND  THE  RESULTANT  TORSOR  OF  THE  FOLLOWING  LOADS//) 
110  F0RMAT(21X*  6HF0RC  ES»  38 X » 1 OHAPPLIED  AT) 

120  FORMAT  (BX,2HFX,1 3X,2HFY, 13X,2HFZ, 1  4X  ,  1  HX,  1 4X  ,  1  HY, 14X,  1HZ) 

130  FORMAT(6(3X,F12.5) ) 

115  FORMAT(24X,7HCOUPLES ) 

125  FORMAT ( 10X»2HMX» 1 5X,2HMY, 15X,2HMZ) 

135  F0RMAT(3(3X,F 14.7  )  ) 


51  FORMAT (5X»24HTHE  RESULTANT  TORSOR  IS,//) 

50  FORMA T { 5 X , 5HF  =  (,  F14.7,7H)  I  +  (,F14.7,7H)  J  +  (,F14.7,3H)  K) 
55  FORMAT ( 5X , 5HM  =  (,  F14.7,7H)  I  +  (,F14.7,7H)  J  +  (,F14.7,3H)  K) 
134  F0RMAT(5X,21HTHE  BAR  REACTIONS  ARE) 

133  F0RMAT(5X,6F14.7) 

136  FORMAT! 1 H  1  ) 

NP  =  0 

READ  5 ,  NPROB 
520  READ  40 

READ  5,  NBAR 

READ  19,  ( (VB( I , J ) , J=1 ,3) , (B( I , JJ ) , JJ=1 ,3) , 1=1, NBAR) 

CALL  STAOET 
READ  1 0 , NFO , NCO 
IF(NFO)  590,600,590 
590  CONTINUE 

READ  20,  ((FORCE!  I  ,  J  )  , J= 1 , 3 ) , ( CORD  (  I , J J ) , JJ  =  1 , 3) , I  *  1 , NFO ) 

600  CONTINUE 

IF(NCO)  610,620,610 
610  CONTINUE 

READ  30,  ((CCUPL(  I , J ) , J= 1 , 3 ) , I = 1 , NCO ) 

620  CONTINUE 

DO  1  I  =  1,6 


ND  =  NOSOL ( I ) + 1 


1  CONTINUE 
CALL  RE  S 1 
DO  200  I  =  1,3 
DO  200  J  =  1 , N6AR 

200  COMX ( I, J)  =  VB ( J, I ) 

DO  201  I  =  1,NBAR 

COMX (4,1)  =  —  V  B { I  ,2)*B(I,3)+VB(I,3)*B( 1,2) 
COMX ( 5 , I )  =  -VB( I  ,3 J *B( I , 1 )+VB( I,  1  )*B( 1 ,3) 

201  COMX (6,1)  =  - VB  C I  ,  1  )*B( I ,2)  +  VB( I ,2)*B( I, 1 ) 
SUM  =  0. 

DO  202  I  =  1,3 
SUM  =  SUM  +  RFCR(  I  1 

202  COMX  (  I , NB  AR+ 1  )  =  RFOR(I) 

co- 

SUM2  =  0. 

DO  203  I  =  4,6 

SUM2  =  SUM2  +  RCOPLl 1-3) 

203  COMX ( I , NB AR+1 )  =  RCOPLl 1-3) 

IF  ( SUM-SUM2 )  613,613,614 

614  ZR  =  (SUM/3. ) *1 .E-10 
GO  TO  615 

613  ZR  =  (SUM2/3. ) el.  E-10 

615  CONTINUE 

IF  { KPROB ( 4 ) )  302,301,302 
301  NUNC  =  6 


CALL  REACT ( COMX , B RE AC , NUNC, Z R , K PROB I  3 1 ) 
GO  TO  2 
302  CONTINUE 

EP  =  .00000001 
IF  INFO)  609,604,609 
609  DO  601  I  =  1,3 
AR  {  I )  =  V  B  C 1  ,  I  ) 

601  BR ( I )  =  VB ( 2 ,  I  ) 

CALL  CROST  (  AR, BR .CROSS  1 
SUM  =  0. 

DO  602  I  =  1,3 

602  SUM  =  SUM  +  RFORI  I ) « RFOR I  I  1 
DEN  =  SQRTF ( SUM ) 

DO  603  I  =  1,3 

603  AR  (  I  )  =  RFORI  D/DEN 
C  =  0. 

DO  606  I  =  1,3 
606  C  =  C  +  CROSS! I ) *  AR I  I ) 

IF  (ABSF(C)-EP)  604,604,605 
605  KPROB (5 )  =1 
GO  TO  2 

604  CONTINUE 

IF  I NCO )  612,611,612 
612  SUM  =  0. 


DO  607  I  =  1,3 

607  SUM  =  SUM  +  RCOPL ( I ) *RCOPL(I ) 

DEN  =  SQRTF(SUM) 

DO  608  I  =  1,3 

608  AR ( I )  =  RCOPL ( I ) / DEN 
CALL  CROST  ( AR , CR OS S , BR ) 

IF( (ABSF(BR( 1 ) )+ABSF( BR(2) )+ABSF( BR{3) ) )-EP)  61 1 ,6 1 1 ,605 
611  CONTINUE 
11  =  1 
J  1  =2 

K1  =3 
NUNC  =3 

■P-  303  DC  304  I  =  1,4 

o 

ps  c  l,  i )  =  coMxm  ,  i ) 

PS ( 2 ,  I)  =  COMXIJl  ,1  ) 

3014  PS  (  3,  I  )  =  C0MXIK1  ,  I  ) 

CALL  REAC1 (PS, B RE  AC, NUNC, ZR,KM) 

IF  (KM)  306,2,306 


3  06 

IF 

(Kl-6) 

308, 

309 

,309 

308 

Kl 

=  K 1  +  1 

GO 

TO  303 

309 

IF 

( Jl-5) 

310, 

311 

,31  1 

310 

Jl- 

-Jl  +  1 

Kl  =  Jl+1 


GO  TO  303 

311  IF  (11-4)  312,313,313 

312  II  =  11+1 

J1  =11+1 
K 1  =  Jl+1 

GO  TO  303 

313  KPROB ( 3 )  =  1 
2  CONTINUE 

PRINT  90 
PRINT  40 
PRINT  101 
PRINT  111 
£  PRINT  121 

PRINT  129, ( <VB( I, J) , J=1 , 3) , (B( I , JJ) , JJ=1 ,3) , 1=1 , NO AR) 
PRINT  91 

IF  ( KPROB ( 1 ) )  157,158,157 

157  PRINT  516 

158  IF  ( KPROB ( 2 ) )  159,165,159 

159  PRINT  517 

165  I F ( KPROB ( 3 ) )  161,  162,161 

161  PRINT  515 

162  IF  ( KPROB ( 4 ) )  163,168,163 

163  PRINT  518 
GO  TO  164 


168  PRINT  521 

164  I F ( KPROB ( 5 ) )  167,166,167 

167  PRINT  91 
PRINT  519 
GO  TO  320 
166  CONTINUE 
PRINT  91 
11=0 

DO  102  I  =  1,6 
II  =  11+1 

IF  { NOSOL ( I ) )  112,102,112 

102  CONTINUE 

-P-  GO  TO  321 

fo 

112  GO  TO  (500,501,502,503,504,5051,11 

500  PRINT  510 
PRINT  91 
GO  TO  320 

501  PRINT  511 
PRINT  91 
GO  TO  320 

502  PRINT  512 
PRINT  91 
GO  TO  320 

503  PRINT  513 


PRINT  91 
GO  TO  320 

504  PRINT  514 
PRINT  91 
GO  TO  320 

505  PRINT  522 
PRINT  91 
GO  TO  320 

321  CONTINUE 
PRINT  100 

IF  (NFC)  316,317,316 

316  CONTINUE 

-P-  PRINT  110 

Vo 

PRINT  120 

PRINT  130,  (  (  FORCE  (I,J),J=1,3),  (CORDl  I ,  J  J  )  ,  J  J=  1 , 3  ),  1  =  1  ,NFO) 

PRINT  91 

DO  560  1=1 , NFO 

DO  560  J= 1 , 3 

FORCE ( I , J ) =0 . 

560  CORD ( I , J ) =0 . 

317  CONTINUE 

IF  ( NCO )  314,315,314 
314  CONTINUE 
PRINT  115 


PRINT  125 

PRINT  135, ( ( COUPL ( I , J J J) , JJJ= 1 , 3) , 1=1 ,NCO) 
PRINT  91 
DO  570  1=1, NCO 
DO  570  J=  1 , 3 

570  COUPL ( I , J  )  =0 . 

315  CONTINUE 

PRINT  51 

PRINT  50, ( RFOR ( J) , J= 1 ,3) 

PRINT  55, (RCOPL(J) , J=l,3) 

PRINT  91 
PRINT  13L 

PRINT  133,  (BREACI  I  )  ,  1=1 , 6) 

320  CONTINUE 
PRINT  90 
DO  16  1=1,5 

16  KPROB ( I )  =  0 
DO  17  I  =  1,6 

17  NOSOL ( I )  =0 

DO  571  I  =  1 , NBAR 
00  571  J  =  1,3 
B( I , J  )  =  0. 

571  VB ( I , J)  =  0. 

DO  572  I  = 


1,6 


PRINT  136 

I F ( NP-NPROB )  520,580,580 


UJ 

ZC 


O si 

N* 

m 


h-  a 

2:  a  a 
01-2: 

O  m  uj 

o 

CO 

m 
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PROGRAM  TESTIS 

DIMENSION  VB(6»3)  »  B ( 6  »  3 ) , KPROB ( 5 ) , NOSOL { 6 ) , RFOR ( 3) ,RCOPL(3), 

1BREAC (6) , COMX  (  6 ,7 ) , AR { 3) , BR{ 3 ) , CROSS ( 3 ) , VP( 5 , 3 ) , PO INK  3, 4 ) , 

2PLA ( 1 2 »  4 ) »  PL ( 4 ) ,PM{4) ,BFORCE(6) ,PS(6,7) 

COMMON  VB  t  B  » NBAR» KPROB, NOSOL , COMX , VP, BFORCE 
5  FORMAT { 15) 

19  FORMAT(3F10.8,3F1  4.7) 

20  FORMAT ( 6F 1 2. 6 ) 

90  FORMAT (/////) 

400 FORMAT ( 80H 

1  ///) 

1 C 1  FORMAT (  5X ,  48HEXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE//) 

■P-  1  1  10F0RMATI 1 1X,25HDIRECTI0N  COSINES  OF  BARS,21X, 

o 

121HPOINTS  OF  APPLICATION) 

121  FORMAT  (8X, 1HL, 12 X, 1 HM, 12X, 1HN, 1 3X, 1HX, 16X, 1HY, 1 6X, 1HZ ) 

129  FORMAT (3(2X, FI  1.7) , 3 ( 3X , F 1 4 . 7 ) ) 

130  F0RMATI5X, 19HALL0WABLE  BAR  LOADS//) 

131  FORMA  T ( 6F 1 4 . 7 ) 

91  FORMAT!//) 

5 1 00F0RMAT { 5X ,  61 HSTRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  C 
1 ONCURRENT ) 

5110F0RMAT15X, 5 9H STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  PA 
1RALLEL) 

512  FORMAT(5X,50HSTRUCTURE  UNSTABLE  INSUFFICIENT  NUMBER  OF  SUPPORTS) 


513  F0RMAT15X » 4 9H STRUCTURE  INDETERMINATE  EXCESS  NUMBER  OF  SUPPORTS) 
S140F0RMAT(5X»63H STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERS 
1 ECT  A  LINE) 

515  FORMAT(5X,27HREACTION  MATRIX  IS  SINGULAR) 

516  FORMAT (5X,28HSUPPORT  POINTS  ARE  ON  A  LINE) 

517  FORMAT l 5X ,  29HSUPP ORT  POINTS  ARE  IN  A  PLANE) 

518  FORMAT (5X, 14HPLANAR  PROBLEM) 

521  FORMAT(5X » 2 OH NOT  A  PLANAR  PROBLEM) 

522  FORMAT! 5X, 51HSTRUCTURE  UNSTABLE  SUPPORTS  ARE  IN  PARALLEL  PLANES) 
61  FORMAT(5X,28HTHE  ALLOWABLE  LOAO  TORSOR  IS) 

50  FORMAT ( 5X  t5HF  =  (,  F14.7,7H)  I  +  (,F14.7,7H)  J  +  (,F14.7,3H)  K) 

55  FORMAT ( 5X » 5HM  =  (,  F14.7,7H)  I  +  (,F14.7,7H)  J  +  (,F14.7,3H)  K) 

•P-  1  36  FORMAT  (  1  H  1  ) 

NP  =  0 

READ  5 »  NPROB 
520  READ  40 

READ  5,  NBAR 

READ  19,  (  (  VB(  I  ,J  ) . J=1 ,3) , ( B( I , JJ) , JJ=1 , 3), 1  =  1, NBAR) 

CALL  STADET 

READ  20,  (BFORCEl  I  )  , I  =  1 , 6 ) 

DO  1  I  =  1,6 
ND  =  NOSOU  I  )  +  1 
GO  TO  ( 1 ,2) ,ND 
1  CONTINUE 


DO  200  J  =  l.NBAR 

200  COMX ( I, J)  =  VB( J, I ) 

DO  201  I  =  1 , NB  AR 

COMX  (  4,1)  =  -VB(I ,2)»B(I,3)+VB( I,3)*B(I,2) 

COMX ( 5 » I )  =  -VB(I  ,3)*B( I, 1 )+VB( I, 1)*B( 1,3) 

201  COMX (6*1)  =  -VB ( I , 1 ) *8 ( I , 2 ) +VB ( I , 2 ) *B( I » 1 ) 

DO  602  I  =  1,3 

RFOR ( I )  =  0. 

602  RCOPL (I)  =0. 

DO  601  I  =  1,3 
DO  601  J  =  1,6 

RFOR ( I )  =  RFOR( I) +BFORCEI J)*COMX( I, J) 

601  RCOPL ( I )  =  RCOPL(  I)+8FCRCE(J)*C0MX( 1+3, J) 

2  CONTINUE 
PRINT  90 
PRINT  40 
PRINT  101 
PRINT  111 
PRINT  121 

PRINT  129, ( ( VB( I , J) , J=1 ,3 ) , ( B ( I , JJ ) , JJ=1 ,3) , 1=1 ,NBAR) 
PRINT  91 

IF  (KPROB(l))  157,158,157 
157  PRINT  516 


158  IF  ( KPRQB ( 2 ) )  159,165,159 

159  PRINT  517 

165  1 F { KP  ROB (3)1  161,162,161 

161  PRINT  515 

162  IF  ( KPROB ( 4 ) )  163,164,163 

163  PRINT  518 
GO  TO  322 

164  PRINT  521 
322  CONTINUE 

PRINT  91 
11=0 

DO  IOC  I  =  1,6 

ii  =  n  +  i 

IF  ( NOSOL (Ill  110,100,110 
100  CONTINUE 
GO  TO  321 

110  GO  TO  (500,501,502,503,504,5051,11 

500  PRINT  510 
GO  TO  320 

501  PRINT  51 1 
GO  TO  320 

502  PRINT  512 
GO  TO  320 

503  PRINT  513 


GO  TO  320 

504  PRINT  514 
GO  TO  320 

505  PRINT  522 
GO  TO  320 

321  CONTINUE 
PRINT  91 
PRINT  130 

PRINT  131 , (BFORCE ( I ) , I  =  1,6) 
PRINT  91 
PRINT  61 

PRINT  50, ( RFOR ( J) , J=  1 ,3) 

PRINT  55, ( RCOPL ( J ) ,J=1 ,3) 

320  CONTINUE 
PRINT  91 
DO  16  1=1,5 

16  KPROB ( I )  =0 
DO  17  I  =  1,6 

17  NOSOL ( I )  =  0 

DO  571  I  =  1 , NB AR 
DO  571  J  =  1,3 
B( I , J)  =  0. 

571  VB( I, J)  =  0. 

NP=NP+ 1 


I F ( NP-NPROB )  520,580,580 
580  CONTINUE 


a 

o  c 

00  UJ 
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APPENDIX  D 


Subroutine  CR03T 


The  subroutine  evaluates  the  cross  product  of  two 
vectors  expressed  in  rectangular  coordinates.  A  Fortran 
expression  is  written  for  each  term  of  the  expansion  of 
CROSS. 


CROSS  s 


A  i  k 


A*  Ay  Az 
Bx  By  Bz 


=  (AyBz-AzBy)i  +  (AzBx-AxBz)j2  +  (AxBy-AyBx)k 
=  CR0SS(l)i  +  CR0SS(2)j  +  CR0SS(3)k 
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SUBROUTINE  CROST  ( AR , BR , CROS S ) 
DIMENSION  AR{3),3R(3),CROSS(3) 
CROSS(I)  =  AR{2)*BR(3)-AR(3) *BR ( 2  ) 
CROSS ( 2 )  =  AR { 3 ) *  BR { 1 ) -AR  C 1 ) *BR ( 3 ) 
CROSS ( 3 )  =  AR( 1 )*BR(2)-AR(2)*BR{ 1 ) 
RETURN 
END 


v_n 

VjO 


APPENDIX  E 

Subroutine  RES1 


The  subroutine  finds  the  resultant  of  a  system  of 
forces  and/or  couples  by  first  adding  up  the  rectangular 
components  of  either  and  then  calculating  the  effect,  if 
any,  of  the  forces  on  the  resultant  couple. 


F  = 

s 

F  • 

etc 

X 

1 

XI 

Mx  = 

£ 

i 

F  . 

XI 

etc 

r  =  rxi  +  ryj_  +  rzk 


F  =  Fxi  +  Fyj  +  Fzk 

M  =  r  x  F 

=  M^i  +  Myj  +  Mzk 


Fxj_  is  the  x  component  of  Ft 
Mx^  is  the  x  component  of  PL 


r  is  the  radius  vector  to  the 
point  of  application  of  the 
force 


rx,  ry  and  ra  are  the  rectan¬ 
gular  components  of  r 

F  is  the  vector  force 

Fx,  Fy  and  Fz  are  the  rectan¬ 
gular  components  of  F 

M  is  the  vector  moment 


M^.,  My  and  Mz  are  the  rectan¬ 
gular  components  of  M 
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SUBROUTINE  RESl 

DIMENSION  V  B  C  6  *  3 )  ,  B ( 6, 3 ) , FORCE! 100,3) ,CORD( 100,3),  C0UPL(50,3)  , 

1 KP  ROB (5) , NOSOL (6)  , RFOR { 3 ) , RCOPL ( 3 ) , BREAC ( 6 ) , COMX { 6 , 7 ) , PS( 6,7 ) , 

2  AR { 3 ) , BR ( 3) , CROSS ( 3 ) , VP ( 5 , 3 ) , PO I  NT ( 3, 4 ) , PL A  ( 1 2 , 4 ) , PL ( 4 ) , PM ( 4 ) 
COMMON  VB,B, FORCE , CORD ,COUPL , NFO, NCO, NBAR , KPROB , NO  SOL , RFOR , RCOPL , 
"l BREAC, COMX, VP 
DO  5  1=1,3 
RCOPL ( I )=  0. 

5  RFOR ( I )=0. 

IF  (NFO)  40,50,40 
40  CONTINUE 
DO  10  1=1,3 
DO  10  J= 1 , NFO 

10  RFOR(I)  =  RFOR ( I ) +  FORCE(J,I) 

50  CONTINUE 

IF(NCO)  60,70,60 
60  CONTINUE 
DO  20  J= 1 , 3 
DO  20  1  =  1  , NCO 

20  RCOPL t J)=RCOPL( J) +COUPLI I, J) 


70  CONTINUE 

IF(NFO)  80,90,80 
80  CONTINUE 

DO  30  1=1 , NFO 

RCOPL ( 1 ) =RCOPL ( 1) -FORCE! 1,2) *CORD { I , 3 ) +FORCE ( I 
RCOPL(2)=RCOPL(2)-FORCE( I,3)*C0RD! I, 1 )+FORCE( I 
30  RCOPL(3)=RCOPL(3)-FORCE! I, 1 )*CORD(I ,2)+F0RCE( I 
90  CONTINUE 
RETURN 
END 


,  3)*CGRD(  1,2) 
,  1  )  *  CORD { 1,3) 
, 2  )*  CORD ( 1,1) 


vn 

O 


APPENDIX  F 


Subroutine  REAC1 

This  subroutine  is  a  modification  of  the  GAUSS2  sub- 

p 

routine  used  by  C.  B.  Bailey  in  his  program  for  solving 
simultaneous  linear  equations.  Whereas  Bailey  was  solving 
the  matrix  equation,  Ax  =  B,  with  a  maximum  of  50  vectors 
B  for  each  coefficient  matrix  A,  nowhere  in  this  program 
is  more  than  one  B  vector  involved  with  a  given  A. 

The  routine  locates  the  largest  first  column  coef¬ 
ficient  and,  if  necessary,  exchanges  the  row  containing 
this  coefficient  with  the  first  so  the  largest  element  is 
A-q  .  Multiples  of  row  one  are  subtracted  from  the  other 
rows  to  make  the  column  one  coefficients  below  the  first 
zero.  Next  the  column  two  coefficients  below  the  first 
are  examined  to  find  the  largest  and  rows  are  exchanged, 
if  necessary,  as  before.  This  Gaussian  elimination  process 
is  continued  until  all  the  elements  below  the  main  diagonal 
are  zero. 

At  each  reduction  step,  the  value  of  the  diagonal 
element  is  compared  with  zero.  A  very  small  diagonal 
element  indicates  ill-conditioned  equations  and  causes 
an  error  output.  This  is  the  "matrix  singularity"  test 
used  in  C0NCURRENT1.  If  the  matrix  is  non-singular,  a 
back-  solving  process  is  used  to  find  the  unknowns. 

2C.  B.  Bailey,  F2  UTEX  LINEQN,  CO  OP  Manual,  Dec.,  I960. 
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SUBROUTINE  REAC1 ( AA, X, N, FP,KERJ 
DIMENSION  A  A  (  6 »  7 ) »X(6) 


KER  =  0 
DO  1»  I  =  l.N 
4  X  (  I  )  =  0. 

NPM  =  N  +  1 
DO  34  L  =  1  ,N 
KP  =  C 
Z  =  0.0 
DO  12  K  =  L.N 

IF (Z-ADSF( AA(K,L) ) )  11,12,12 

11  Z  =  A  B  S  F  (  A  A  (  K  ,  l ) ) 

KP  =  K 

12  CONTINUE 
IF(L-KP)  13,20,20 

13  DO  14  J  =  L , NPM 
Z  =  A  A  (  L  ,  J  ) 

AA(L,J)  =  AA( KP, J ) 

14  AAIKP.JJ  =  Z 

20  IF{ ABSF(AA(L,L) )-FP)  50,50,30 

30  IF  (L-N)  31,40,40 

31  LP1  =  L+1 

DO  34  K  -  LP1 , N 

IF  ( A A ( K , L ) )  32,34,32 


32  RATIO  *  AA(K,L)/AA(L*L) 

DO  33  J  =  LP1.NPM 

33  AA (  K ,  J  )  =  AA<K,J)  -  RAT I 0* AA { L , J ) 

34  CONTINUE 

40  DO  43  I  =  1  i N 
II  =  N  +  1  -  I 
S  =  0.0 

IF  (I  I  —  N  )  4  1,43,43 

41  I  I P 1  =  II  +  1 

DO  42  K  =  I  I  P 1 , N 

42  S  =  S  +  AA( II , K ) *  X ( K ) 

*1 3  X  (I  I  )  =  (  AM  I  I  ,  NPM)  -S  )  /  AA  (II  ,  I  I  ) 
RETURN 
50  KER  =  1 
RETURN 


APPENDIX  G 


Subroutine  PLANE 

Given  a  vector  with  direction  cosines  1,  m  and  n 
and  one  point,  the  problem  is  to  find  a  pair  of  planes 
whose  intersection  is  parallel  to  the  vector  and  passes 
through  the  point . 

First,  the  subroutine  checks  to  see  if  the  vector 
is  "nearly"  parallel  to  a  coordinate  axis.  "Nearly"  is 
defined  as  having  a  direction  cosine  greater  than  0.#, 
i.e.,  lying  within  about  37°  of  the  coordinate  directions. 
If  the  vector  is  "nearly"  parallel  to  one  axis,  it  is 
crossed  with  vectors  along  the  perpendicular  axes. 

If  the  vector  is  not  "nearly"  parallel,  a  check  is 
made  to  see  whether  it  lies  within  about  6°  of  a  coordinat 
plane,  i.e.,  the  direction  cosine  is  less  than  0.1.  If 
it  is  close  to  a  coordinate  plane,  one  of  the  crossing 
vectors  is  chosen  perpendicular  to  that  plane. 

If  neither  of  the  above  tests  is  met,  the  given 
vector  is  crossed  with  vectors  in  the  x  and  z  directions. 

There  are  now  two  cross  product  vectors  of  the  form: 
X^  ~  AjL  +  B_j  +  Ck 

Since  these  vectors  are  perpendicular  to  the  planes  con¬ 
taining  the  vectors  which  formed  the  cross  products,  the 
dot  product  of  X^  anc*  an  arbitrary  vector  in  one  of  the 
planes  is  zero.  An  arbitrary  vector  is  of  the  form: 
(x-x0)i  +  (y-y0)j.  +  (z-z0)k 
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where  (x0,yo,zo)  is  the  known  point.  The  dot  product  has 
the  form: 

A(x~xo)  +  B(y~yo)  +  C(z-zo)  =  0 
This  may  be  rearranged  to: 

Ax  +  By  +  Cz  =  D 
where:  D  =  Ax0  +  By0  +  Cz0 

The  coefficients  A,  B,  C,  and  D  are  stored  in  the  array 
PLA. 
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SUBROUTINE  PLANE! CC , DO , PL , PM ) 

DIMENSION  CC ( 3 ) ,QO(3),PL(4),PM{4),AA(3),BB(3) 


DO  7  K  =  1,3 
AA ( K )  =  0  . 

PL ( K )  =  C. 

-  PM(K)  =  0. 

7  P.B(K)  =  0. 

PL ( 4 )  =  C. 

PM  ( *i )  =  C. 

DO  1  K  =  1  , 3 
IFCAOSFCCC(K) )-.B)  1,1,3 

1  CONTINUE 

DO  1 1  K  =  1 ,3 
IF ( ABSF (CC ( K) )-. 1 )  9,9,11 
11  CONTINUE 

5  AA (  1  )  =  1  . 

B  B  (  3  )  =  1. 

GO  TO  4 

6  AA ( 1 )  =  1  . 

BB ( 2 )  =  1. 

GO  TO  4 

2  A  A  (  2  )  =  1  . 

BB ( 3  J  =  1  . 


4  CONTINUE 

CALL  CROSUCC,  AA,  PL) 

CALL  CROST(CC,BB»  P.M  ) 

8  DO  10  K  =  1,3 

PL (4)  =  P  L ( 4 )  +  PL(K)*DD(K) 
10  P M ( 4  >  =  P M ( 4 )  +  PM ( K  )  * 00 ( K ) 

RETURN 

3  GO  TO  (2,5,6) ,K 

9  GO  TO  (6,2,5) ,K 
END 


o 


APPENDIX  H 


Problem  Format 


TEST1A 


Fortran  Format 

Symbol(s) 

15 

NPROB 

Hollerith 

15 

NBAR 

3F10.S,3F14.7 

VB ,  B 

13,12 

NFO ,NC0 

6F12.6 

FORCE, CORD 

3F14.7 

COUPL 

Input  data 

Total  number  of  problems 
Problem  title  -  up  to  SO 
characters 
Number  of  links 
Direction  cosines,  coordi¬ 
nates  of  links 
Number  of  forces,  couples 
Force  components,  points  of 
application 
Couple  components 


Repeats  from  the  Hollerith  title  for  each  problem* 
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TEST1B 


Fortran  Format 

Symbol(s) 

Input  Data 

15 

NPROB 

Total  number  of  problems 

Hollerith 

Problem  title  -  up  to  BO 

characters 

15 

NBAR 

Number  of  links 

3F10.S,3F14.7 

VB ,  B 

Direction  cosines,  coordi¬ 

nates  of  links 

6F12.6 

BFORCE 

Allowable  link  loads 

Repeats  from  the 

Hollerith 

title  for  each  problem 
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APPENDIX  I 


STADET  Flags 


KPROB  =  Kind  of  Problem 

1  Support  points  are  on  a  line 

2  Support  points  are  in  a  plane 

3  Reaction  matrix  is  singular 

4  Planar  problem  (if  4  is  1) 

Not  a  planar  problem  (if  4  is  0) 

5  Resultant  load  is  not  in  plane  of  structure 


N0S0L 

1 

2 

3 

4 

5 

6 


No  solution  because: 

Structure  unstable  -  lines  of  action  of  supports 
are  concurrent 

Structure  unstable  -  lines  of  action  of  supports 
are  parallel 

Structure  unstable  -  insufficient  number  of 
supports 

Structure  indeterminate  -  excess  number  of  supports 
Structure  unstable  -  lines  of  action  of  supports 
intersect  a  line 

Structure  unstable  -  supports  are  in  parallel 
planes 


66 


APPENDIX  J 


STADET  List 
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o 

ca 


SUBROUTINE  STADET 

DIMENSION  VB  (  6 , 3)  , B ( 6 , 3 ) , FORC E ( 1 00 , 3 ) , CORD ( 1 00, 3 1 ,  C0UPL(50,3)  , 

1 KPROB ( 5 ) , NOSOL (6) ,RF0R(3) ,RC0PL(3) , BREAC ( 6 ) , COMXt 6 , 7 ) , PS ( 6, 7 ) , 

2  AR ( 3 ) , BR ( 3  J .CROSS (3) , VP ( 5 , 3 ) , PO I  NT ( 3, 4 ) , P L A  (1  2 , 4 ) , PL ( 4 ) , PM ( 4 1 
COMMON  VB,B, FORCE , CORD, COUPL, NFO.NCO, NBAR , KPROB , NO  SOL , RFOR, RCOPL , 
1 BRE  AC , COMX, VP 
C  COP  L ANAR  1 

SUM  =  0. 

DO  401  I  =  1,3 
DO  401  J  =  1,3 
4  01  SUM  =  SUM  +  B( I , J  ) 

EP  =  ( SUM/9. ) * ! • E-6 
NBA  =  NBAR  -1 
DO  1 1  I  =  1 , NBA 
DO  11  J=  1 , 3 

11  VP(I,J)=B((I+1),J)-B(1,J) 

I  =  1 
J  =  0 
Jl  =  0 

.  23  IF  (ABSF ( VP( I , 1 ) ) +ABSF(VP( I ,2) )+ABSF( VP( I ,3) )-EP)  22,22,21 

22  1=1+1 

IF  ( I -NBAR  )  23,315,23 
21  IF  (I-NBAR+1)  24,25,25 
24  K  =  I  +  1 


26 

32 

33 

25 

41 

42 

43 

31 

5 

51 

61 

152 

156 


IF  (ABSFIVPIK,  1  ))+ABSF(VPIK,2) ) +ABSF ( VP ( K , 3 ) )-EP)  32,32,31 
I  F  (  K  +-  1-NBAR  )  33,2  5,  25 
K=K+  1 
GO  TO  26 

IF  (NEAR-3)  41,42,43 

NOSOL (3)=1 

RETURN 

KPROB I  1  )  =  1 

KPR03I4)  =  1 

GO  TO  74 

NO SOL  <4)  =  1 

RETURN 

DO  5  L  =  1,3 
ARIL)  =  V  P (  I, L) 

BRIL)  =  VPIK.L) 

CALL  CROST  ( A R  ,  BR , CROSS ) 

IF  (ABSFICROSSl 1)  ) + ABS F I C ROS S I  2 ) ) + ABS F I C ROS S I  3 ) )-EP)  51,51,152 
IF  IK  +  l-NBAR)  61,25,25 
K  =  K  + 1 
GO  TO  26 

I F I NBAR-3 ) 4 1 , 156,  52 
KPRGBI2)  =  1 
KPROBI4)  =  1 
GO  TO  73 


52  J  =  J  +  l 

IF  ( J- I  )  53,52,53 

53  IF  ( J-K)54,52,54 

54  C  =CRCSS ( 1 ) *VP ( J , 1 )+CR0SS(2 )*VP(J,2)+CR0SS( 3  )  *  VP  ( J ,3) 
J1  =  J1  +  1 

IF  (ABSF(C)-EP)  62,62,81 
62  IF  ( J 1 +  3-NBAR  )  52,72,72 

72  KPR08 {  2  )  =  1 

1  =  1 

73  C=CROSS (  1  )*VB( 1,1  ) +CROSS { 2 ) * VB (  1,2) +CROSS ( 3 ) * V 8 ( 1,3) 
IF  (ABSF(C)-EP)  32,82,31 

81  IF  CJCAR-6)  83,76,05 

85  NO  SOL ( 4 )  =  1 
RETURN 

83  NOSCl (  3  )  =  1 
RETURN 

82  IF  (  I-NBAR  )  06,87,  87 

86  1=1+1 

GO  TO  73 

87  I F ( NB AR-3  )  41,08,43 

88  KPROB ( 4 )  =  1 


C  PL ANAR  1 

74  1=2 

79  DO  6  L  =  1,3 
AR ( L )  =  V  B (  1  ,L) 

6  BR( LJ  =  VB ( I , L ) 

CALL  CROST  (  AR,  83.  »  CROSS  ) 

101  IF  (ABSF(CROSSM) ) + ABSF ( CROSS ( 2 ) ) +ABSF ( CROSS ( 3) ) -E P )  103,103,102 

103  I  F (  1-3  )  104,105,1  05 

104  1=3 

GO  TO  79 

105  NOSOL (2 )  =  1 
RETURN 

102  IF  ( 1-3)  106, 176, 176 

106  C  =  CROSS ( 1  ) *VB ( 3 , 1  )  +  CROSS ( 2 ) * VB I  3 , 2 ) +CROSS ( 3 ) * VB ( 3, 3 ) 

IF  (ABSF(C)-EP)  77,77,41 

176  C=CROSS ( 1  ) *VB ( 2, 1  ) +CROSS ( 2 ) *VB( 2, 2 ) +CROSS ( 3 ) * VB ( 2,  3) 

IF  (ABSF(C)-EP)  77,77,41 


C  PARALLEL1 

76  1  =  1 
109  1=1+1 
K  =  0 
M  =  0 

DO  7  L  =  1 , 3 
ARIL)  =  V  B ( 1 y  L ) 

7  BR ( L )  =  VB ( I, L) 

CALL  CROST  ( AR, RR , CROSS ) 

IF ( ABSF ( CROSS ( 1 ) )  + ABS F ( CROSS ( 2 )  )  +  ABS F ( CRO S S ( 3 )  ) - EP  )  103,108,501 

108  IF  (  I-.NBAR)  109,11  1,111 
111  NOSOL (2) =1 
RETURN 

501  C  =  0. 

K  =  K+  1 

DO  502  L  =  1,3 
C  =  C  +  CROSS ( L 1  *  VB ( K , L ) 

502  CONTINUE 

IF  (ABSF(C)-EP)  503,503,504 


IF  (K-NBAR)  501,505,505 
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C  PARALLLL2 

38  1  =  1 
J  =  2 

1 1  =0 

116  DO  8  L  =  1,3 

AR ( L )  =  VB(  I , L  ) 

8  8R(L)  =  VB(J,L) 

CALL  CROST  ( AR , BR , CRC SS ) 

I F ( ABSF (CROSS (11) + ABSF (CROSS (2 ) ) + ABSF ( CROSS ( 3 ) ) - EP )  1 13, 1 13, 1 12 
113  11=11+1 

112  I F  (  J-6)  1  1  >4,  115,1  15 

1 1 U  J=J+1 

GO  TO  life 

115  I F (  1-5)  118,117,117 
118  1=1+1 
J  =  1  +  1 

GO  TO  116 

117  I F (  11-6)  77,315,315 


CONCURRENT! 

77  LL  =  0 
402  LL  =  LL  ♦  1 
DO  1  N  =  1,3 
AR(N)  =  VB(LL,N) 

1  B  R  (  N  )  =  D ( LL , N ) 

CALL  PLANE  ( A R , BR , P L , PM ) 

DO  316  I  =  1,4 
KO  =  2*  LL 
K I  =  KC— 1 

PLA  ( K I , I )  =  PL ( I ) 

316  PLAIKC.I)  =  Pf'.(I) 

IF  (LL-N6AR)  402,2,2 

2  CONTINUE 
K 1  =0 
K2  =  4 

204  KL  =  4 

DO  3  K  =  1,4 
DO  3  L  =  1  ,4 

3  COMX ( K, L  )  =  PLAIK  +  K1 ,L) 

CALL  RE AC  1 (COMX,BREAC,KL,EP,KM) 
IF  (KM)  201  ,202,201 

202  IF ( K2-2*NBAR)  203,205,205 

203  K 1  =  Kl+2 


K2  =  K2  +  2 
GO  TO  204 

205  RETURN 
201  CONTINUE 

KL  =  3 

11  =  1 

12  =  2 
13  =  3 

208  00  206  L  =  1,4 

PS(  1,11  =  COMX ( 1 1  , L ) 

PS ( 2  ,  L )  =  COMX I  12 ,L) 

PS ( 3 , L )  =  COMX (13 , L 1 

206  CONTINUE 

CALL  R E AC  1 ( PS , BRE AC , KL , EP , KM ) 
IF(KM)  207,4,207 

207  CONTINUE 
13  =  13+1 

IF  (13-4)  208,208,209 

209  12  =  12+1 
13  =  12+1 

IF  (13-4)  208,208,210 

2  10  II  =  I  1 +1 
12  =11+1 
13  =  12+1 


208,208,202 


IF  (13-4) 

4  CONTINUE 
15=0 
16  =  0 

213  15  =  15+1 
DO  =  0. 

DO  314  K  =  1,3 

314  DO  =  DO  *  BREAC(K )»PLA( 15, K) 
RID  =  PLA (15,4)  -  DD 

IF  ( ABSF(RID)-FP)  312,312,311 

312  16  =  16+1 

311  I F ( 1 5-2  BAR )  213,317,317 

317  IF  (  NBAR-3  )  318,318,313 

318  IF  (16-6)  202,321,321 
321  NOSOL ( 1 )  =  1 

RETURN 

313  IF  (16-8)  202,315,315 

315  NOSOL ( 5 )  =  1 
RETURN 

END 


APPENDIX  K 


Problem  Types 

In  this  appendix  "bars"  and  "supports"  are  used  as 
synonyms  for  "links".  See  Section  3»  Testing  of  Program, 
for  a  discussion  of  the  problems. 
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1.  UNSTABLE  STRUCTURE  -  COINCIDENT  POINTS  OF  SUPPORT 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

-.7388637 

.6562276 

.1531198 

16.5000000 

1 1.7000000 

3 . 2500000 

.7388648 

.6562286 

.1531200 

16.5000000 

i  11.7000000 

3.2500000 

-.7388637 

.6562276 

.1531198 

16.5000000 

11.7000000 

3.2500C00 

NOT  A  PLANAR  PROBLEM 

STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT  A  LINE 


-o 

vO 


2.  PLANAR  PROBLEM 


TOO  FEW  SUPPORTS 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


Z 

6.  400C0CO 
1  7.  630CCC0 


NOT  A  PLANAR  PROBLEM 


DIRECTION  COSINES  OF  BARS 
L  M  N 

,0591533  .4050570  .6543264 

0915723  .9030155  .4197378 


POINTS  OF  APPLICATION 
X  Y 

1.3000000  25.2000000 

3.7500000  1.0400000 


STRUCTURE  UNSTABLE  INSUFFICIENT  NUMBER  OF  SUPPORTS 


3. REDUNDANT  STRUCTURE 


MORE  THAN  THREE  BARS 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

3 F  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

2 

1  . ooooocc 

.0000000 

.0000000 

.0000000 

. OOOOOCO 

. oooccco 

.oooooco 

1  .coooooc 

. 0000000 

.oooooco 

.oococco 

. oooccco 

. oooooco 

.0000000 

1 .0000000 

.0000000 

.0000000 

.oooooco 

1.0000000 

.0000000 

.ocooooo 

10.0000000 

.0000000 

. oooooco 

. ocooooo 

1  .cocoooo 

. 0000000 

10.0000000 

. 0000000 

.OOOOCC  0 

.0000000 

.0000000 

1.0000000 

10. ooooocc 

.OOOOOCC 

.oooccco 

NOT  A  PLANAR  PROBLEM 

STRUCTURE  INDETERMINATE  EXCESS  NUMBER  OF  SUPPORTS 


A.  UNSTABLE  THREE  BAR  STRUCTURE 


BARS  ARE  PARALLEL 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

.oooooco 

1 .OOCOOOO 

.OOOOOOO 

.ooooooc 

.OOOOOCO 

. oooooco 

. OOOOOCO 

1 .OOCOOOO 

.0000000 

10.0000000 

.0000000 

. oooooco 

.OOOOOCO 

-1 .OOCOOOO 

.0000000 

10.0000000 

.0000000 

. OOOOOCO 

SUPPORT  POINTS  ARE  ON  A  LINE 
PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LIMES  OF  ACTION  OF  SUPPORTS  ARE  PARALLEL 


5.  UNSTABLE  THREE  BAR  STRUCTURE 


BARS  ARE  NOT  IN  A  PLANE 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


01 

L 

. 1 ooocco 

. 9000CC0 
.4246410 


RECTICN  COSINE 
M 

.9000000 
.  100C000 
.1000000 


S  OF 


BARS 

N 

424264 1 
424264  1 
9000000 


POI 

X 

8.9000C0C 
0 . 9000C00 
16.5000000 


OF  APPLICATION 
Y 

8.90CCCCO 
8 . 9C00CC0 
1  1 .  1000000 


NTS 


z 

5. OCOCCCO 
5.00000C0 
10. 10000C0 


SUPPORT  POINTS  ARE  ON  A  LINE 
PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  INSUFFICIENT  NUMBER  OF  SUPPORTS 


6.  UNSTABLE  THREE  BAR  STRUCTURE 


LINES  OF  ACTION  OF  SUPPORTS  INTERSECT 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

1.0000000 

.oocoooo 

.0000000 

1 .2000000 

3.6000000 

7.2000000 

.0000000 

1  .0000000 

.0000000 

1  .2000000 

3.6000000 

7.2000000 

.4472136 

.8944272 

.0000000 

3.2000000 

7.6000000 

7.2000000 

SUPPORT  POINTS  ARE  ON  A  LINE  O 

PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  CONCURRENT 


OQ. 

■F- 


7.  REDUNDANT  STRUCTURE 


MORE  THAN  THREE  BARS 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTI CN  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

1 .ooocooo 

.oooooco 

.0000000 

.OOCOOOO 

. CCOOOOO 

. ocoocco 

1 .0000000 

.0000000 

.OOOOOCO 

lO.OCOOOCO 

.oooooco 

. oooocco 

1 .0000000 

.ooocooo 

.0000000 

10.0000000 

.oooooco 

10.0000CCC 

1 .oooooco 

.oocoooo 

.oooooco 

20.0000000 

.oooocco 

10.0C0CCCC 

1 .coococo 

.ooccoco 

.0000000 

10.0C0000C 

.cccooco 

-10. oooccco 

1 .oooooco 

.oocoooo 

. ocooooo 

20. OOCOOOO 

.0000000 

-10. oooooco 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 


STRUCTURE  INDETERMINATE  EXCESS  NUMBER  OF  SUPPORTS 


OQ. 

vn 


8.  UNSTABLE  STRUCTURE 


SUPPORTS  ARE  PARALLEL 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L  K 

N 

X 

Y 

z 

-l.OOOOOCO  .OOOOOOO 

. OOOOOOO 

10.1000000 

8.9000000 

5. OOOOCCO 

1.0000000  .0000000 

.0000000 

8.9000000 

8 . 9000CC0 

5.000CCCC 

1.0000000  .oooooco 

.OOOOOOO 

1 6.5000000 

1  1  .  1000000 

5. OOOOCCO 

SUPPORT  POINTS  ARE  IN  A 
PLANAR  PROBLEM 

PLANE 

STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  PARALLEL 


9.  UNSTABLE  STRUCTURE  -  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

.OOOOOOO 

1.0000000 

.0000000 

1.2000000 

3.6000000 

7.2000CC0 

.4472136 

.8944272 

.0000000 

3.2000000 

7.6000000 

7 . 2000CC0 

1.0000000 

.0000000 

.0000000 

3.6000000 

3.6000000 

7.2000000 

SUPPORT  POINTS  ARE  IN  A  PLANE 
PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  ARE  CONCURRENT 


oa 

-^3 


10.  UNSTABLE  STRUCTURE  -  LESS  THAN  SIX  BARS 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


D  I 

RECTICN  COSINES 

OF 

BARS 

POINTS 

OF  APPLICATION 

L 

N 

N 

X 

Y 

z 

1.0000000 

.0000000 

. 0000000 

6. 3000000 

8.300CCCC 

. cooocco 

.0000000 

.0000000 

1 

.0000000 

8. 7000CCC 

7. 30CCCC0 

.oooccco 

. 0000000 

.0000000 

1 

.ccooooo 

5 . 2  OCOOOO 

8. 5CC0C00 

. OOOOCCO 

. 0000000 

1 .coccooo 

.ccooooo 

9.3000000 

2.  1000000 

. OOOOCCO 

.0000000 

1 .0000000 

.0000000 

5.4000000 

3. 20000C0 

. oooccco 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 

STRUCTURE  UNSTABLE  INSUFFICIENT  NUMBER  OF  SUPPORTS 


oa 

oa 


11.  REDUNDANT  STRUCTURE 


MORE  THAN  SIX  BARS 


EXAMINE  THE 

stability  of 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POI  NTS 

CF  APPLICATION 

L 

M 

N 

X 

Y 

Z 

-1 .0000000 

.oocooco 

1.0000000 

9. 100000C 

5.4000CCC 

3 . 20000C0 

1  .oooocco 

.0000000 

. 0000000 

6.30000C0 

8. 30000C0 

. OOOOCLC 

1.0000000 

.0000000 

.OCOOOOO 

8.7000C00 

7. 3000C00 

. OOOOCCO 

-1 .0000000 

.0000000 

.0000000 

9. 1000000 

6.4000CC0 

.OOOOCCO 

-1 .0000000 

.cocoooo 

.0000000 

5.20000CC 

8 . 500CCCC 

. OOOOCCO 

. 0000000 

1  .cocoooo 

. OCOOOOO 

9 . 30000C0 

2.  1000000 

. oooocco 

.0000000 

1  .0000000 

.OCOOOOO 

5.4000000 

3.2000000 

. oooocco 

NOT  A  PLANAR 

PRCCLEM 

STRUCTURE  INDETERMINATE  EXCESS  NUMBER  OF  SUPPORTS 


12.  UNSTABLE  STRUCTURE 


SUPPORTS  ARE  PARALLEL 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

V 

N 

X 

Y 

z 

.0000000 

1 

.0000000 

.0000000 

.0000000 

.0000000 

. oooocco 

.0000000 

-1 

.ooooocc 

.0000000 

10.0000000 

.oooccco 

.OCCGCCO 

.0000000 

1 

.oocooco 

.0000000 

.0000000 

.0000000 

1 c. oooocco 

.0000000 

-1 

.0000000 

.0000000 

10.0000000 

. 0000000 

10. oooocco 

.0000000 

1 

.0000000 

.coococo 

.0000000 

. oooooco 

-10. oooccco 

.0000000 

-1 

.cocoooc 

.0000000 

10.0000000 

. oooocco 

-1C.00CCCC0 

SUPPORT  POINTS  ARE  IN  A 

PLANE 

NOT  A  PLANAR 

PROBLEM 

STRUCTURE 

UNSTABLE  LINES 

OF  ACTION  OF 

SUPPORTS  ARE  PARALLc 

L 

vO 

O 


13.  UNSTABLE  STRUCTURE 


SUPPORTS  ARE  IN  PARALLEL  PLANES 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

1.0000000 

.0000000 

.0000000 

.0000000 

.0000000 

.oooooco 

.0000000 

1 .0000000 

.0000000 

.oooooco 

.oooooco 

. oooccco 

1 . oooooco 

.0000000 

.0000000 

10.0000000 

.oooooco 

10.CC0CCC c 

.oooccco 

1.0000000 

.0000000 

10.00000CC 

. oooooco 

1 0. oooccco 

.0000000 

1 .coooooo 

.0000000 

10.0000000 

.0000000 

-  1 0. oooooco 

1.0000000 

.0000000 

.0000000 

10.0000000 

. 0000000 

-10. oooooco 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  SUPPORTS  ARE  IN  PARALLEL  PLANES 


14.  UNSTABLE  STRUCTURE 


FOUR  OF  SIX  SUPPORTS  ARE  PARALLEL 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF 

BARS 

PCINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

l 

-1 . oooooco 

.0000000 

.0000000 

9.1 OCOOOO 

6.4000000 

. OOOOCC  0 

.0000000 

1 .0000000 

.0000000 

6.3000000 

8. 3000000 

. oooocco 

. 0000000 

.0000000 

1 

.0000000 

8.7000000 

7.3000000 

.oooooco 

.0000000 

.oocoooo 

1 

.0000000 

5.20000CG 

8 . 50CC000 

. oooocco 

. 0000000 

.0000000 

-1 

.0000000 

9.3000000 

2. 1000000 

. oooocco 

.0000000 

.0000000 

1 

.0000000 

5 . 4000000 

3.2000000 

. oooocco 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT  A 


LINE 


15.  UNSTABLE  STRUCTURE 


LINES  OF  ACTION  OF  SUPPORTS  INTERSECT 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

K 

N 

X 

Y 

z 

.4472136 

.8944272 

.OOOOOOO 

3.2000000 

7.6000000 

7. 200CCC0 

.OOOOOOO 

.8944272 

-.4472136 

1.2000000 

7.6000000 

5.200CCC0 

-.4472136 

.8944272 

.OOOOOOO 

-.8000000 

7.6000C00 

7. 200CCC0 

.0000000 

.8944272 

.4472136 

1  .2000000 

7.6000000 

9. 20000C0 

I . OOOOOOO 

.OOOOOOO 

.OOOOOOO 

1.2000000 

7.6000000 

7. 2000CC0 

-1.0000000 

.OOOOOOO 

.OOOOOOO 

1.2000000 

7.6000000 

7.200CCC0 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT  A  LINE 


16.  UNSTABLE  STRUCTURE 


LESS  THAN  SIX  BARS 


EXAMI NE 

THE  STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DI 

RECTIGN  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

V 

N 

X 

Y 

z 

1.0000000 

.0000000 

.0000000 

8.7000000 

7. 3000000 

. OOOOCLC 

-.5801556 

-.1)850556 

.6543245 

1 . 30000CC 

25. 2000CC0 

6.4CCCCCC 

.0915722 

-.90301 45 

.4197373 

3.75C00C0 

1 .0400000 

1  7. 630C0C0 

. 9000000 

.L2U2641 

-.1 OOOOCO 

1  1 . 65C0000 

12.9500000 

1  6.  300CCL0 

.9000000 

.  1000000 

.424264  1 

. OOCOOCC 

.OOOOCCO 

.OOOOCCO 

NOT  A  PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  INSUFFICIENT  NUMBER  OF  SUPPORTS 


17.  REDUNDANT  STRUCTURE  -  MORE  THAN  SIX  BARS 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

2 

.  4242641 

.9000000 

. 1000000 

.0000000 

16.5000000 

1 1 . 100CCC0 

.0915722 

-.9030145 

.4197373 

3.750000C 

1 . 0400CC0 

1  7 . 630CCC  0 

-.5801556 

-.4850556 

.6543245 

1 . 3000C00 

25.2COOOCO 

6.4000CC0 

1.0000000 

.0000000 

. OCOOOOO 

.0000000 

.OCOOOOO 

. OOOOOCC 

1 .0000000 

.0000000 

.0000000 

6. 300.0000 

8. 3000CC0 

. OOOCCLO 

-1 .0000000 

.0000000 

.0000000 

9. 1000CCC 

6. 40CCCC0 

. OOOCCLO 

. 9000000 

.1000000 

.424264  1 

16.50C00C0 

1  1  .  1000000 

10.  1000CC0 

NOT  A  PLANAR  PROBLEM 


STRUCTURE  INDETERMINATE  EXCESS  NUMBER  OF  SUPPORTS 


18.  UNSTABLE  STRUCTURE 


SUPPORTS  ARE  PARALLEL 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECT 

ION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

V. 

N 

X 

Y 

z 

. 0000000 

1 .0000000 

.0000000 

.0000000 

4. T  OOOOCO 

. oooocco 

.0000000 

1 .0000000 

.0000000 

TO. OOOOOCO 

. OOCCCCO 

.  oooccco 

.oooocco 

1 .0000000 

.0000000 

10.0000000 

4 . 1 000000 

1C.  oooocco 

. 0000000 

1 .0000000 

.oooooco 

20.0000000 

. OCOOOOO 

1 0. oooocco 

. oooooco 

1 .0000000 

.0000000 

1 0. oooooco 

. OOOOOCO 

-1 0. oooooco 

.0000000 

1 .0000000 

.coooooo 

20.0000000 

.OOOOCCO 

-10. OOCCCCO 

NOT  A  PLANAR 

PROBLEM 

STRUCTURE  UNSTABLE  LINES 

OF  ACTION  OF 

SUPPORTS  ARE  PARALLE 

L 

19.  UNSTABLE  STRUCTURE 


SUPPORTS  ARE  IN  PARALLEL  PLANES 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECT 

ION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

1 . oooooco 

.0000000 

. 0000000 

1 .6500000 

2. 9500000 

6. 3000CC0 

1 . ooococo 

.0000000 

. ococooo 

12.0000000 

29.0000000 

. oooccco 

-  1 . coooooo 

.00-00000 

.0000000 

12.0000000 

29 . OOOCCCC 

. oooccco 

.  1 000  oco 

.90C0CC0 

.424264  1 

.OOOOCCC 

. OCCCCCO 

. oooccco 

-  1 . oooooco 

.cocoooo 

. 0000000 

8. 75C0000 

. ccoocco 

6. 66000CO 

. 1000000 

.9000000 

. 4  2U  26A  1 

8.7500000 

. oooooco 

6. 6600CC0 

NOT  A  PLANAR 

PROBLEM 

STRUCTURE  UNSTABLE  SUPPORTS  ARE  IN  PARALLEL  PLANES 


20.  UNSTABLE  STRUCTURE  -  FOUR  OF  SIX  SUPPORTS  ARE  PARALLEL 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECT 

ION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

Z 

-1 .0000000 

.0000000 

.OOOOOCO 

.0000000 

.OOOOOCO 

2 

1 . 300CCC0 

1 . oooooco 

.0000000 

.0000000 

.0000000 

.OOOOOCO 

2 

1 . 3000CC0 

-1 .0000000 

.0000000 

.OOOOOCO 

5.2000000 

8. 5C00000 

. OOOOOCO 

.0000000 

1 .0000000 

.0000000 

9. 300000C 

2. 100CCC0 

. OOOOOCO 

.ocooocc 

.0000000 

1 .0000000 

5. 4000CCC 

3. 200CCC0 

. cooccco 

1.0000000 

.0000000 

.0000000 

.0000000 

.OOOOOCO 

. oooooco 

NOT  A  PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT  A  LINE 


21.  UNSTABLE  STRUCTURE  -  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECT 

ION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

1.0000000 

.0000000 

.0000000 

1.2000000 

3.6000000 

7.20000C0 

.0000000 

1 .0000000 

.0000000 

1.2000000 

3.6000000 

7. 2000CC0 

.4472136 

.8944272 

.0000000 

3.2000000 

7.6000000 

7.200CCC0 

.0000000 

.8944272 

-.4472136 

1.2000000 

7. 6000CC0 

5.200CCCO 

-.4472136 

.8944272 

.0000000 

-.8000000 

7.6000000 

7 . 2000CC0 

.0000000 

.8944272 

.4472136 

1.2000000 

7.6000000 

9. 2000CC0 

NOT  A  PLANAR  PROBLEM 


STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS  INTERSECT  A  LINE 
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22a.  PLANAR  PROBLEM 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


D I 

RECTI CN  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

N 

N 

X 

Y 

2 

1  .0000000 

.coccooc 

.0000000 

21  .2000000 

.0000000 

. oooocco 

I  . 0000000 

.0000000 

.0000000 

35.3000000 

21.2000000 

. oooocco 

.0000000 

-1 .0000000 

.0000000 

35.3000000 

21 .2000000 

. oooocco 

SUPPORT 

POINTS  ARE  ON  A 

L  INE 

PLANAR  PROBLEM 


FIND  THE  RESULTANT  TORSOR  OF  THE  FOLLOWING  LOADS 


FX 

FORCES 

FY 

FZ 

X 

APPLIED  AT 

Y 

Z 

-692. 17000 

737.21000 

.00000 

16. 10000 

6.78000 

.occco 

17.02000 

7092.22000 

.00000 

9.70000 

5.95000 

.ooccc 

1021 . 12000 

223.24000 

.00000 

3.40000 

3.920C0 

.oocco 

COUPLES 

MX 

MY 

MZ 

. 0000000 

.0000000 

2  59 . OOCOOOO 

. 0000000 

.0000000 

52.4900000 

.0000000 

.0000000 

739.5000000 

THE 

RESULTANT  TORSOR 

IS, 

F  = 

(  345.9700C00) 

I  +  (  8052. 

6699998)  J  +  ( 

.OOCOOOO ) 

K 

M  = 

(  .0000000) 

I  +  ( 

OOOCOOO)  J  +  ( 

77956.8741 989) 

K 

THE 

BAR  REACTIONS  ARE 

-9385.2741888  9731. 

2441885  -8052.6699998 

.OCGOCOO 

•CCCGOCO 

.occoccc 

22b.  PLANAR  STRUCTURE  WITH  IMPROPER  LCAD 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POINTS 

OF  APPLICATION 

L 

N 

N 

X 

Y 

Z 

.  lOOOOCO 

.9000000 

.4242641 

1 .3500000 

7. 6500000 

19. 730CCC  0 

.  lOOOOOO 

. 90CC0C0 

.4242641 

1 . 3500CCC 

7 . 650CCCC 

1 9. 730CCC0 

. 1000000 

.9000000 

.424264  1 

9 . 7500CCC 

7.6500000 

19. 7300CC0 

SUPPORT  POINTS  ARE  ON  A  LINE 
PLANAR  PROBLEM 


RESULTANT  LOAC  IS  NOT  IN  PLANE  OF  STRUCTURE 


102 


23a.  OBLIQUE  PLANE  PROBLEM  RESULTANT  FORCE  IS  IN  PLANE 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L  M 

N 

X 

Y 

z 

. 7308637  -.6562276 

-.1531198 

1 . 3000000 

25.2000000 

6 . 400CCC0 

-.5001556  -.4350556 

. 6543245 

16.5000000 

1 1 .7000000 

3. 2500CC0 

-.0915722  .9030145 

-.4197373 

3.7500000 

1 . 0400000 

17. 6300CC0 

SUPPORT  POINTS  ARE  IN  A 
PLANAR  PROBLEM 

P  LANE 

FIND  THE  RESULTANT  TORSO*  OF  THE  FOLLOWING  LOADS 


FORCES 

FX  FY 

-73.88637  65.62276 

-58.01556  -45.50556 

9.15722  -90.30145 

FZ 

15.31198 

65.43245 

41.97373 

X 

1 . 300C0 
16.5CCC0 

3. 75000 

APPLIED  AT 

Y 

25. 200C0 

1  1.70000 

1 . 04000 

Z 

6. 4 COCO 
3.250C0 
17. 630C0 

THE  RESULTANT  TCRSOR  IS, 

F  =  (  -  122. 7447C60)  I  +  ( 

M  =  (  2524.7461075)  I  +  (  - 

-73.1842410)  J  + 
1756.9239447)  J  + 

(  122. 71  81  56C  )  K 

(  1477.5325581)  K 

THE  BAR  REACTIONS  ARE 

-100.0000000  100.0000000 

-ICO. 0000000 

.0000000 

.OOOCCCO 

.OOOOCCC 

23b.  OBLIQUE  PLANE  PROBLEM  RESULTANT  FORCE  IS  NOT  IN  PLANE 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIPECTI 

ON  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

7388640 

.6562279 

. 1 531 199 

1 6. 5000CCC 

1  1  . 7C0CCC0 

3.250CCC0 

5801557 

.L3SC556 

.6543246 

3 .7500000 

1 . 0400000 

17.6300CL0 

0915722 

.9020145 

.4197373 

1.3000000 

25.2000000 

6.  4000CL0 

SUPPORT  POINTS  ARE  IN  A  PLANE 
PLANAR  PROBLEM 


RESULTANT  LOAD  IS  NOT  IN  PLANE  OF  STRUCTURE 


23c.  OBLIQUE  PLANE  PROBLEM  RESULTANT  FORCE  IS  NOT  EXACTLY  IN  PLANE 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

.7383637 

-.6562276 

-. 1 531 1 98 

1 . 30000CC 

25.20CCCC0 

6.4000CC  0 

-.5801556 

-.4850556 

.6543245 

1 6. 50000C0 

1  1.7000000 

3. 2500000 

-.0915722 

.9030145 

-.4 197373 

3.75C0000 

1 .0400000 

17. 6300CCO 

SUPPORT  POINTS  ARE  IN  A  PLANE 
PLANAR  PROBLEM 


RESULTANT  LOAD  IS  NOT  IN  PLANE  OF  STRUCTURE 


H 

O 

-p- 


24.  STABLE  SPATIAL  PROBLEM 


EXAMINE  THE 

STABILITY  OF 

THE  FOLLOWING 

STRUCTURE 

DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

z 

.  1000000 

. 9000000 

.  424264  1 

.0000000 

.oooocco 

.oooooco 

.9000000 

. lOCOOOO 

.4242641 

.0000000 

.0000000 

. oooooco 

.9000000 

.4242641 

-.  1000000 

1  1.6500000 

12.9500000 

16. 300CCC0 

. 1 000000 

.9000000 

-.4242641 

11.6500000 

12. 9500000 

1 6.300CCC0 

.4242641 

. 10C0C00 

.9000000 

.ooooccc 

9.7000CC0 

7 . 900CCC  0 

.4242641 

.9000000 

. 1 COOOOO 

.0000000 

9.7000000 

7.9000GC0 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 


FIND  THE  RESULTANT  TORSOR  OF  THE  FOLLOWING  LOADS 


FX 

FORCES 

FY 

FZ 

X 

APPLIED  AT 

Y 

Z 

84.27500 

-10. 10000 

19.24000 

-6.5C000 

5.  14000 

1 .50CL0 

COUPLES 

MX  MV  MZ 

78 . 3C00C00  68.5000000  87.1100000 


THE  RESULTANT  TORSOR  IS, 


F  =  (  84.2750000)  I  +  (  -10.1000000)  J  +  ( 

M  =  (  192.3436000)  I  +  (  319.9725000)  J  +  ( 


THE  BAR  REACTIONS  ARE 

39.0014903  82.3859488  44.5109861  -28 


19.2400000)  K 
-280.4135000)  K 


6600889  -39.1484822 


-47.3504682 
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25.  STABLE  S P AT  I AL  . PRCBL EM 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 
L  M 

OF  BARS 

N 

X 

POINTS  OF  APPLICATION 

Y 

Z 

. 9000000 

. 1000000 

.4242641 

.0000000 

1  1.6500000 

.oooocco 

.  1 000000 

.9000000 

.4242641 

1  1 .6500000 

.CCOOCCO 

.cooccco 

.  1000000 

.9000000 

-.4242641  - 

21 .1200000 

1  .5000000 

.oooocco 

. 9000000 

.4242641 

-.1000000 

.0000000 

.0000000 

16. 30000C0 

. 4242641 

. 1000000 

.9000000 

.0000000 

9. 70000C0 

7. 900CCCC 

.424264 1 

.9000000 

. 1000000 

7. 50000CC 

9. 7000CC0 

.OOOOCCO 

NOT  A  PLANAR  PROBLEM 

FIND  THE  RESULTANT  TORSOR  OF  THE  FOLLOWING  LOADS 


FORCES 

APPLIED  AT 

FX 

FY 

FZ 

X 

Y 

30.  12000 

39.79000 

1 0.20000 

-8.70000 

2. 13000 

MX 

19.9000000 


COUPLE  S 
MY 

90C. 1500000 


MZ 

10.2900000 


THE  RESULTANT  TORSOR  IS, 


F  =  (  30.1200000) 

M  =  (  1.0360000) 


I  +  (  39. 790C000 )  J 

I  +  (  1019.0100000)  J 


+  (  10.20C0000)  K 

+  (  -400.0386000)  K 


Z 

l.OOCCO 


THE  BAR  REACTIONS  ARE 

-121.0548607  -106.1962450 


167. 0425823 


-32.8508693 


41.9365735 


89.0738227 
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26.  PLANAR  PROBLEM 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTI 

CN  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

Z 

1 

.0000000 

.0000000 

. 0000000 

21.2000000 

.oooocco 

.OOOCCCC 

1 

.0000000 

.0000000 

.0000000 

35.30CCCC0 

21 . 2CCCCC0 

.cccccto 

.OCOCOOO  -1 

.ooccooo 

.0000000 

35.3000000 

2  1 . 2C00C00 

. oooocco 

SUPPORT  POINTS  ARE  ON  A  LINE 
PLANAR  PROBLEM 


ALLOWABLE  BAR  LOADS 


-9385.2741873  9731.2441878-8052.6699986  .OOCOCCC  .OOGCCCC  .CCCCCCC 


THE  ALLOWABLE  LOAD  TCRSOR  IS 

F  =  (  345.9700003)  I  +  (  8052.6699986)  J  +  (  .OCCOCCO)  K 

M  =  (  .0000000)  I  +  (  .0000000)  J  +  (  77956.8741760)  K 


£0T 


27.  OBLIQUE  PLANE  PROBLEM  RESULTANT  FORCE  IS  IN  PLANE 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTION  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

Z 

.7388637 

-.6562276 

-.1 531198 

1 .3000000 

25.2000000 

6.400CCC0 

— .5801556 

-.4850556 

.6543245 

16.5000000 

1  1 .7000000 

3. 2500CC0 

-.0915722 

.9030145 

-.4  197  37  3 

3 . 75C0000 

1.C400CC0 

1  7.630CCC  0 

SUPPORT  POINTS  ARE  IN  A 
PLANAR  PROBLEM 

PLANE 

ALLOWABLE 

BAR  LOADS 

-lOC.COOCOOO 

1C0.C000000 

-1 00.0000000 

.0000000 

. 0000000 

.CCOOOOO 

THE  ALLOWABLE  LOAD  TORSOR  IS 

F  =  (  -  122.7  447C60 )  I  +  (  -73.1842410)  J  +  ( 
M  =  (  2524. 7461 C74 )  I  +  (  - 1 756. 9239446 )  J  +  ( 


122.7181 560 )  K 
1477.5325581  )  K 
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STABLE  SPATIAL  PROBLEM 


EXAMINE  THE  STABILITY  CF  THE  FOLLOWING  STRUCTURE 


DIRECTI 

ON  COSINES 

OF  BARS 

POINTS  OF  APPLICATION 

L 

M 

N 

X 

Y 

Z 

.  1000000 

.9000000 

.4242641 

.0000000 

.0000000 

.oooocco 

. 90000C0 

.1000000 

.4242641 

.0000000 

. cooooco 

.ococccc 

.9000000 

.424264 1 

-. 1000000 

1  1 . 65C0CCC 

1 2 . 9500CC0 

16.3C0CCC0 

.  10CC000 

.9000000 

-.424264  1 

1  1 . 65CCCC0 

12.9500000 

16. 300CCC0 

.4242641 

.  10C0000 

. 9CC00C0 

.0000000 

9.7000000 

7.9000CC0 

.4242641 

.9000000 

.  1C00000 

.0000000 

9. 70CCCC0 

7.900CCCC 

SUPPORT  POINTS  ARE  IN  A  PLANE 
NOT  A  PLANAR  PROBLEM 


O  ALLOWABLE  BAR  LOADS 


39.0014903  82.3859488  44.5109861  -28.6600889  -39.1484822  -47.3504682 


THE  ALLOWABLE  LOAD  TORSOR  IS 

F  =  (  84.2750000)  I  +  {  -10.100CC00)  J  +  (  19.24COCOO)  K 

M  =  (  192.3436002)  I  +  {  319.9724995)  J  +  (  -28C . 4  1 34 998 )  K 


no 


29.  STABLE  SPATIAL  PROBLEM 


EXAMINE  THE  STABILITY  OF  THE  FOLLOWING  STRUCTURE 


DIRECTI 

L 

CN  COSINES 
M 

OF  BARS 

N 

X 

PCINTS  OF  APPLICATION 

Y 

Z 

. 9000000 

.1000000 

.4242641 

.0000000 

1  1 .6500000 

.OOOOCCO 

. 1000000 

.9000000 

.424264 1 

1  1 .6500000 

.OOOOCCO 

. oooccc  c 

. 1000000 

.9000000 

-.424264 1 

21. 120C00C 

1 . 500CCCO 

.ooocccc 

.9000000 

.4242641 

-.1000000 

.OCCCCOC 

. CCCCCCO 

16. 3000CC0 

.4242641 

. 1 COOOOO 

. 9000000 

.0000000 

9.7C00CC0 

7.9000CC0 

.4242641 

.9000000 

.  1COOOOO 

7.5000000 

9. 7000000 

. OOOOCCO 

NOT  A  PLANAR  PROBLEM 


ALLOWABLE  eAR  LOADS 


-121 . C548607  -1C6. 1962450  -32.8508693  41.9365735  89.0738227  167.C425823 


THE  ALLOWABLE  LOAC  TORSOR  IS 

F  =  (  30 . 1 20CC00 )  I  +  {  39.7900000)  J  +  (  1O.20CO0CC)  K 

M  =  (  1.8360000)  I  +  (  1019.0100006)  J  +  (  -400.0385996)  K 


APPENDIX  L 


Definition  of  ’’Zero” 

The  computer  must  convert  any  non-integer  number  to 
an  approximate  binary  form  which  will  fit  its  word  length. 
The  conversion  and  length  limitation  allow  an  accuracy  of 
about  ten  decimal  digits.  With  each  calculation  round-off 
errors  accumulate.  Whereas  the  flow  charts  ask  if  a 
quantity  is  zero,  what  must  really  be  considered  is 
whether  the  quantity  is  less  than  the  maximum  expected 
round-off  error.  If  there  are  no  conflicting  requirements, 
the  maximum  round-off  error  may  be  estimated  and,  with 
perhaps  some  leeway,  defined  as  "zero".  However,  where 
different  orders  of  magnitude  may  be  involved  in  the  same 
problem  or  in  subsequent  problems  or  where  it  is  advisable 
to  have  a  large  tolerance  for  one  calculation  and  a  small 
tolerance  for  another,  the  situation  becomes  far  more 
complicated.  One  cannot  define  the  "zeros"  in  terms  of 
the  number  of  significant  figures  in  an  expression;  a 
specific  magnitude  is  required. 

An  attempt  was  made  to  tie  the  size  of  the  "zero"  to 
the  size  of  the  structure.  It  was  felt  that  barring  un¬ 
reasonable  complications,  such  as  having  a  very  small 
structure  at  a  great  distance  from  the  origin^  the  system 
of  averaging  the  coordinates  of  the  points  of  application 
and  multiplying  this  average  by  the  computer  accuracy 
would  give  a  suitable  zero.  However,  instead  of  being 
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able  to  use  10  ^  times  the  average  coordinate,  it  was 
necessary  to  use  a  "zero"  10,000  times  as  big  to  get 
proper  answers.  Further  investigation  revealed  that 
the  difficulty  probably  involved  the  use  of  subroutine 
REAC1  in  C0NCURRENT1 . 

In  subroutine  REAC1  a  small  value  of  "zero"  allows 
small  diagonal  elements  which,  during  the  back-solving 
process,  generate  large  coordinates  of  the  point  of 
intersection.  When  these  large  values  are  then  used 
in  C0NCURRENT1,  the  round-off  error  becomes  much  greater 
than  the  "zero"  defined  above.  Unsuccessful  efforts  were 
made  to  modify  the  "zero"  obtained  above  in  a  logical  way 

for  use  in  C0NCURRENT1  and  REAC1 . 

-6 

A  multiplier  of  10~  gives  proper  answers,  at  least 
for  this  set  of  problems,  but  there  is  no  other  justifi¬ 
cation  for  its  use.  The  author,  regretfully,  did  not 
have  sufficient  time  to  pursue  this  investigation. 
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